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ON FORCE-FREE MAGNETIC FIELDS 


By 8S. CHANDRASEKHAR 
UNIVERSITY OF CHICAGO 
Communicated November 10, 1955 


1. Introduction —Liist and Schliiter' have recently pointed out that cosmic 
magnetic fields might often satisfy the condition 


curl H X H = 0, (1) 


where H denotes the intensity of the magnetic field. The reason for this, essen- 
tially, is that on account of the high electrical conductivity of stellar material of 
even low density, large currents can flow in the space surrounding the stars; the 
usual assumption (appropriate only for a vacuum) that the current density j 
should vanish outside the photospheres of stars ceases to be even an approximate 
description of the true state of affairs. Under the circumstances existing, it would 
be more accurate to say that the prevailing magnetic field is such that the Lorentz 
force 


l 
£=jxH=7 ulHxH (2) 


vanishes. The presence of the magnetic field will not, then, affect the motions, 
and a stationary state can be envisaged. 

Fields which satisfy equation (1) will be called force-free. 

The force-free condition (1) clearly requires that , 


cul H = oH, (3) 


where a is some function of position. The case when H has symmetry about an 
axis and a is a constant is of considerable practical interest. It is the case which 
Liist and Schliiter have studied in some detail; but they have failed to notice 
that for this case the problem admits an explicit solution. It is the object of this 
paper to present this solution. 

2. The Basie Differential Equations~-A magnetic field none of whose com- 
ponents are azimuth-dependent can always be expressed as the sum of a toroidal 
and a poloidal field in the form 


H = 1, X rT + curl (1, X rP), (4) 


where 1, is a unit vector in the direction of the axis of symmetry (which we shall 
assume to be the z-axis) and P and 7 are any two functions which do not depend 


1 


| 
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on the azimuthal angle g. [Note that for H defined as in equation (4) the condi- 
tion div H = 0 is identically satisfied.]| We shall consider P and 7 as functions 
of either r and & (in spherical polar co-ordinates) or @ and z (in cylindrical polar 
co-ordinates). 

In cylindrical polar co-ordinates, H given by equation (4) takes the form 


oP 10 sp 
H= let ofl, + 1, (5) 


where lg and 1, are unit vectors in the radial (@) and the transverse (¢g) directions. 
For H given by equation (5) it can be readily verified that 


oT 
curl H = —o ly — oA;P1, + 1., (6) 


(7) 
is the operator defining the axially symmetric wave equation in five dimensions. 
From a comparison of equations (5) and (6) it follows that condition (3) requires 

a grad oP = grad oT (8) 
and 
AP = (9) 


A further condition is that grad w*P and grad a are in the same direction. This 
can be deduced as follows: Taking the divergence of equation (3), we obtain 
grad a: H = 0. (10) 


Substituting for H from equation (5) and remembering that grad a has components 
only in the w- and the z-directions, we have 


Oa O Oa O 
(oP) 0. (11) 


An equivalent form of this last equation is 
grad oP X grad a = 0; (12) 


and this is the condition stated. 
3. The Case a = Constant.—In this case equation (12) is clearly fulfilled, while 
equation (8) requires that 


= wT + constant. (13) 


From equation (5) it is evident that an additive constant in #*P has no physical 
consequence. We may therefore assume without loss of generality that 


aP =T. (14) 


Equation (9) now becomes 


2 
where 
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= (15) 


In spherical polar co-ordinates this equation for P is 


40 1 — dp O ) 


where = cos 

Separable solutions of equation (16) exist, and the general solution can be ex- 
pressed in terms of the fundamental solutions 

e rT 3 
p, = my, (17) 

where C,"(u) denotes the Gegenbauer polynomials? and @,, + s/,(ar) is a general 
cylinder function of order n + */2 and represents an arbitrary linear combination 
ot the Bessel functions J, + »,,(ar) and J + s/,)(ar). 

With the solution for the poloidal function P given by equation (17), we can 
write the corresponding magnetic field in the form (cf. eq. [5]) 


H = a(l, X r)P, — — (1. X 7) X grad P,o. (18) 
In spherical polar co-ordinates the various components of H are given by 


4. The Boundary Conditions—The conditions to be satisfied at an interface 
where @ has a jump are (1) that H should be continuous and (2) that the normal 
component of the current density 7 should also be continuous. The latter condi- 
tion requires the continuity of the normal component of curl H. Since curl H = 
aH, the second condition requires the continuity of the normal component of 
aH, On the other hand, the first condition requires that the normal component 
of H be continuous. Since @ changes discontinuously at an interface, the two 
conditions can be simultaneously met only if the normal component of H vanishes. 
The conditions at an interface are, therefore, 


H,=0 and H, and Hy, are continuous. (20) 


From equation (19) it follows that the vanishing of H, at r = R (say) can be 
met only if 


C,, + 0, (21) 


i.c., if aR is a zero of the chosen cylinder function. If the region in question 
extends from r = R, tor = Re, then the linear combination of J, + 1/,and J + 1/4) 
should be such that it vanishes at r = R; andr = Rs. Such cylinder functions for 
various initially assigned values of R:/R, have been explicitly characterized in 
another connection.® 

According to equation (19), the g-component of H vanishes, together with H,, 
on an interface. Finally, the continuity of H, requires that 


3 
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(22) 


be continuous on an interface. Since the foregoing is a linear combination of 
€,, + /,(ar) and a@’, + s/,(ar) and the former vanishes for r = (cf. eq.[{ 21])it is 
clear that the required condition is 


a’, + s/,(ar) is continuous at an interface. (23) 


It can be verified that the conditions found by Liist and Schliiter in their paper 
refer (in our present terminology) to the special case n = 0. 
5. The Energy tn the Magnetic Field—A quantity of interest is the energy in 


the magnetic field, 
1 +1 
| H|%? dr dude 
Jr, J-1 Jo 


1 +1 
| dr dy. (24) 


Since the toroidal and poloidal parts of H are at right angles to each other, 
|H|? = or? sin? P,? + | curl (1, X (25) 


The contribution to F from the toroidal part is readily evaluated. We have 


R2 +1 
E; = ; a? arré*, + sar) du (1 — pw?) (26) 
Ri 
Using the orthogonality property of the Gegenbauer polynomials‘ and of the 
cylinder functions, we have 
1 mn + 1)(n + 2) 
Ep = 4. tp 27 
where Nj, ,, + «/, is the normalization integral for the orthogonal functions provided 
by the functions @,, 4 s/,(a, r) (ef. eq. [3] of an earlier paper’). 
To obtain the contribution to & from the poloidal part, we proceed as follows: 
We have 


1 
Ep = fff eurl(1, rP,,)-eurl (1, rP,,) dx; dirs. (28) 


By an integration by parts we get 
Ep = — f fia. xX curl (1, X rP,)+dS + 
8r Ss 
fff (1, X r)P, curl [eurl (1, rP,,) | dx; daz (29) 
The surface integral vanishes on account of the boundary conditions which require 


P,, to be zero on the bounding surfaces (ef. eq. [21]). The remaining volume 
integral becomes (cf. eqs. [6] and [15]) 
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ip = — 11, X r\2P,AsP,, dx; dz dis 


2 
fff — dr dude. (30) 


It is thus seen that the contribution to the magnetic energy by the toroidal part and the 
contribution by the poloidal part are exactly equal. ‘Thus 
1 1 2 
B = 2B, = 2B, = +a! a + +2) , 
2%(n + */2)(P(?/2) 
6. The Equations for the Lines of Force.—The lines of force are determined by 
the equations 


(31) 


dr rdu — u*) de 


H, H, 


One integral of these equations (namely, that which follows from the first two) 
can be readily found. We have 


(32) 


+ 1Aar)(L — = constant. (33) 


The research reported in this paper has in part been supported by the Geophysics 
Research Directorate of the Air Force Cambridge Research Center, Air Research 
and Development Command, under Contract AF 19(604)-299 with the University 
of Chicago. 

'R. Liist and A, Schliiter, “Kraftfreie Magnetfelder,”’ Z. Astrophys., 34, 263-282, 1954. 

2 The notation is that of G. N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge: 
At the University Press, 1944); see pp. 50 and 369. 

3S. Chandrasekhar and Donna Elbert, “The Roots of J — Ji41/,(dn) X 
J—y+1/(d) = 0," Proc. Cambridge Phil. Soc., 49, 446-448, 1953. 

*Cf. A. Sommerfeld, Partial Differential Equations of Physics (New York: Academic Press, 
1949), p. 232. 


THE EQUILIBRIUM OF MAGNETIC STARS 
By 8S. CHANDRASEKHAR AND Kevin H. PRENDERGAST 
UNIVERSITY OF CHICAGO 
Communicated November 10, 1955 


1. Introduction.—By considering a particularly simple example, Fermi and 
Chandrasekhar! showed that magnetic stars cannot, in general, be spherical. 
Later investigations? carried out with a view to establishing more concrete results 
have proved inconclusive. In this paper we shall prove a theorem which charac- 
terizes in a very general way the magnetic field which can prevail in an incom- 
pressible fluid in equilibrium under its own gravitation. A less ec nplete charac- 
terization for compressible fluids is also given. The discussion is limited to con- 
figurations with an axis of symmetry. 
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2. The Problem.—The equation governing the hydrostatic equilibrium of a 
configuration with a prevalent magnetic field is 

grad p = pgrad V + &, (1) 


where p denotes the pressure, p the density, V the gravitational potential, and £ 
the Lorentz force related to the current density j j 


(2) 


In the usual hydromagnetic approximation, 


1 
j= — 1H 
i= 


£= curl H X H. 
do 
Taking the curl of equation (1), we obtain 
1 
— grad p X grad V = curl £ = Pr curl [curl H X Hi). (5) 
+7 


If we restrict ourselves to configurations with an axial symmetry (about the 2- 
axis), so that p and V are azimuth-independent, the vector on the left-hand side 
has a nonvanishing component only in the transverse ¢g-direction. Consequently, 
curl £ can have no nonvanishing component in the radial (@) and in the z-directions. 
Thus 


(curl £)g = (curl £), = 0. (6) 
If the fluid is incompressible, 
curl £=0, (7) 
and we must in addition require that 
(curl £), = 0. (8) 


The problem now is: What do conditions (6) and (7) imply for the magnetic field? 

3. The Characterization of the Magnetic Field—-We shall suppose that none of 
the components of the prevailing magnetic field is azimuth-dependent. We can 
then express the magnetic field, quite generally, as the sum of a toroidal and a 
poloidal field in the ferm (for the notation see the preceding paper*) 


H = (1, X r)T + curl (1, X rP), (9) 


or, in cylindrical polar co-ordinates, 


oP 1o 
H lg + + (@’P) 1,. 


The curl of this field is given by 


| 
and 
a0) 
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oT 10 
curl H = — lw — oA;P1 
de + (@ ) 2) (11) 
where it may be recalled that A; is the Laplacian operator for an axially symmetric 
function in five-dimensional Euclidean space. The corresponding Lorentz force 
is given by 


4rf A3P OP) ly + 


oro... \ 
7 i+ AsP71., (12) 


and the curl of £ is given by 
oL 


url £ + (= =) 1+ (@£,)1 
cur — — — 


The vanishing of the # and the z-components of curl £, therefore, requires that 


OL, 
Oz 


= 0 and ras (@L£,) = 0. (14) 


These conditions can be met by a nonsingular £, only when 
£, = 0. (15) 
According to equation (12), this implies 


271) — 2P\ — 2 = 0 
( ) 5, (@*P) (@*P) 


_ 
(17) 


We have thus proved that in an axially symmetric configuration in gravitational 
equilibrium the prevalent magnetic field should be such that 


oT = (18) 


where F is a function of the specified argument. 

If we suppose that the configuration is in addition incompressible, then the o- 
component of curl & should also vanish; this further condition gives (ef. eqs. 
{12] and [13]) 

—A,P — (wP) — T— (@&T)¢ —A =0. (19 


On simplifying this equation, we find that we are left with 


(AP, oP) _ 


The general solution of this equation can be found as follows: 


(20) 


7 
or 
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First, we observe that we can define a function G(@?P) such that 
O(a*P, G(@*P)/a*) oT” 
O(a, 2) 


For, on expanding the Jacobian on the left-hand side, we find that the equation 
reduces to 


(21) 


G— = — oT. 
2 5, OP) (22) 
Since #7 is a function of oP only, the foregoing equation is equivalent to 
d 
2G(@*P) = ST"): 2 
(@’P) (23) 


this defines the function G which will satisfy equation (21). 
Returning to equation (20), we can combine it with equation (21) in the following 
manner: 
2) 


(24) 


In other words, 


1 
+ G(@’P) = (25) 


where # is an arbitrary function of the argument. Using equation (23), we can 
rewrite equation (25) in the form 


+ &(@*P), (26) 


d 


which makes the relation between P and T explicit. 
Equation (26) characterizes the most general magnetic field which can prevail in an 
axially symmetric configuration of an incompressible fluid in gravitational equilibrium. 
If 7 = 0, equation (26) reduces to 
= %(@°P); (27) 
this is equivalent to a condition first obtained by Ferraro,‘ who restricted himself 
to poloidal fields. On the other hand, if we let 
T = aP and 0, (28) 


the equation reduces to 
= —a’P; (29) 
this is the equation characterizing force-free fields discussed in the preceding paper. 
The simplest generalization of force-free fields applicable to the problem on hand 
is obtained by setting 
T = aP and © = constant, (30) 
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in which case 
= — + constant. (31) 


The properties of the resulting equilibrium configurations will be investigated 
separately. 

4. An Application to Hydrodynamics.—It is well known that the equations of 
motion of an incompressible inviscid fluid can be written in the form 


—u X curlu = —grad (? + (32) 
ol p 2 


Accordingly, for stationary flows, the velocity field must satisfy the condition 
curl (u X curlu) = 0. (33) 


This is formally the same as the equation for H (cf. eqs. [5] and [7]) considered 
in this paper. Consequently, if we restrict ourselves to axially symmetric station- 
ary flows and represent the velocity field as the sum of a toroidal and a poloidal 
field in the manner of equation (9), then the theorem included in equation (26) 
is applicable to the hydrodynamical problem. These remarks solve generally 
a problem recently considered by Liist and Schliiter.’ 


The research reported in this paper has in part been supported by the Geophysics 
Research Directorate of the Air Force Cambridge Research Center, Air Research 
and Development Command, under Contract AF 19(604)-299 with the University 


of Chicago. 


1. Fermi and 8. Chandrasekhar, “Problems of Gravitational Stability in the Presence of a 
Magnetic Field,’ Astrophys. J., 118, 116-141, 1953. 

2V. C. A, Ferraro, “On the Equilibrium of Magnetic Stars,’ Astrophys. J., 119, 407-412, 
1954; G. Gjellestad, “On Equilibrium Configurations of Oblate Fluid Spheroids with a Magnetic 
Field,” Astrophys. J., 119, 14-33, 1954; P. H. Roberts, “The Equilibrium of Magnetic Stars,” 
Astrophys. J., 122, 508-512, 1955. 

Chandrasekhar, “On Force-free Magnetic Fields,’ these Procrepinas, 42, 1-5, 1956. 

*V.C. A. Ferraro, “On the Equilibrium of Magnetic Stars,”’ Astrophys. J., 119, 407-412, 1954. 

5 R. Liist and A. Schliiter, “Eine spezielle Art nichtwirbelfreier Lésungen der hydrodynamischen 
Gleichungen,” Z. angew. Math. u. Mechanik, 35, 45-47, 1955. 
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CHANGES IN THE COMPOSITION OF AN ACTINOMYCIN COMPLEX 
DURING GROWTH OF A STREPTOMYCES CULTURE 
By W. A. Goss, E. Karz, and 8. A. WAKSMAN 
INSTITUTE OF MICROBIOLOGY, RUTGERS UNIVERSITY, STATE UNIVERSITY OF NEW 
JERSEY, NEW BRUNSWICK 
Communicated November 1, 1955 


The actinomycins thus far isolated from the culture media of certain species of 
Streptomyces have been shown to consist of a group of closely related chemical com- 
pounds composed of a quinonoid chromophore linked to a peptide.’ 2 Physical 
and chemical differences observed among the various actinomycins so far de- 
scribed are considered to be due to variations in the number, arrangement, and 
kinds of amino acids present in the peptide;? the chromophoric moiety is believed 
to be identical in each of these substances.’ The actinomycins reported in the 
literature have been designated as A, B, C, D, I, J, and X.4 Paper chromatography 
and countercurrent distribution of these actinomycins have revealed that they 
are generally composed of a mixture of several components present in different 
proportions.” 4 

Whereas a number of investigations have been carried out on the isolation and 
characterization of the actinomycin complexes and on the classification and identi- 
fication of the microérganisms producing these substances, insufficient attention 
has been given to the influence of nutrition and conditions of cultivation on the 
nature of the complex formed by a given organism. Brockmann observed? that 
the proportion of various components of an actinomycin complex was modified to 
some extent by changes in the nitrogen source of the organism. For example, 
when KNO; was supplied to actinomycin X—producing strains, production of 
actinomycin X; was slight, whereas with glycine the level of actinomycin X, in the 
complex increased considerably. However, in none of Brockmann’s studies did a 
variation of cultural conditions cause a minor component to become a major one. 
Such a reversal in the ratio of the components of an actinomycin complex has 
been found during growth of a Streptomyces culture and is the subject of the present 
communication. 

Streptomyces 3720, isolated from garden soil, initially was found to produce an 
actinomycin complex of the B type during growth on a soya peptone—glucose 
medium. In a subsequent experiment conducted under comparable environ- 
mental conditions, the actinomycin complex formed was identified as the I type. 
The sole difference between the two experiments was the length of incubation of 
the culture. In the first case the material was extracted from 5-day-old cultures 
and in the second from 7-day-old cultures. To resolve this discrepancy, a more 
comprehensive investigation was carried out. 

Experimental.—Five milliliters of a 48-hour-old shake culture of Streptomyces 
3720 was inoculated into 2-liter Erlenmeyer flasks each containing 500 ml. of soya 
peptone medium.® The flasks were then incubated at 28° C. on a rotary shaker. 
At 12-hour intervals aliquots of the broth were taken for the extraction of actino- 
mycin. Using circular paper chromatography (solvent system: 5 per cent sodium 
orthocresotinate and ethyl acetate—n-butyl ether [2:1]), a qualitative determina- 
tion of the complex was achieved. The individual zones on a paper chromatogram 
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were then cut out, the actinomycin eluted with acetone-water (9:1), and the rela- 
tive percentage of each component in a complex calculated after photometric de- 
termination of the eluate with a Beckman DU spectrophotometer. 

It can be seen (Fig. 1) that with an increase in time of incubation there is a con- 
siderable change in the relative percentage of the components found in the actino- 
mycin complex. After 2 days’ incubation the complex was of the B type; there- 
after, a gradual transition took place, and by 5 days the complex was definitely of 
the I type. 


100 


= 
WwW 
2 
= 
re 
2 
a 
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Fic. 1.—Relative per cent of the components in the actinomycin 
complex produced by Streptomyces 3720. The dashed line repre- 
sents the slow-moving component, the dot-dashed line represents 
the intermediate component, and the solid line represents the 
fast-moving component on a circular paper chromatogram using a 
5 per cent sodium orthocresotinate poll ethyl acetate/n-butyl ether 
system. 


Similar studies were performed with Streptomyces 3720, employing a chemically 
defined medium consisting of 1 per cent galactose, 0.2 per cent glutamic acid, 0.1 
per cent K:HPO,, 0.0025 per cent MgSO.*7H.O, 0.0025 per cent FeSO4-7H,0, 
0.0025 per cent ZnSOy-7H,O, and 0.0025 per cent CaCl,-2H,O in distilled water. 
Once again the transition from the B- to the I-type complex occurred. 

To determine whether different strains within the parent culture were responsible 
for the formation of the B- and I-type complexes, several colonial types found in the 
parental culture were isolated and grown separately in soya peptone—glucose 
medium. Analysis of the actinomycin produced during growth of these isolates 
revealed that all initially formed the B-type, and subsequently the I-type, complex. 
Additional studies of the parent culture as well as of the variant types derived from 
it are in progress. On the other hand, a preliminary investigation of Streptomyces 
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3686, an actinomycin I producer, revealed the fact that only the I-type complex 
was formed throughout the entire incubation period. 

Unpublished data’ on the physical and chemical properties of pure components of 
certain actinomycin complexes (B, D, I) have revealed that some of these are identi- 
cal. Therefore, differences among certain actinomycin complexes may depend 
solely on the relative proportions of comparable components. As demonstrated in 
this preliminary report, such proportions can be modified to a considerable extent 
during prolonged cultivation of an organism. What factors may be operating to 
bring about this change are unknown at the present time. It remains to be deter- 
mined whether this modification might involve an effect of the gradually changing 
cultural environment upon the physiological activities of the organism or whether 
such changes may influence the selection of specific types with different synthetic 
abilities.® 

A detailed study of the culture Streptomyces 3720 revealed it to be a typical repre- 
sentative of the S. antibioticus group. It is to be recalled that it was a member of 
this group of actinomycetes that was first found’ to produce actinomycin, later 
designated as “actinomycin A.” 

Summary.—Streptomyces 3720, a member of the S. antibioticus group, produces a 
significant change in the relative percentage of the components present in the 
actinomycin complex during growth on either a complex organic or chemically de- 
fined medium. This organism produces in the medium during the early stage of 
growth the actinomycin B complex. On continued incubation of the culture, the 
actinomycin B is changed to the actinomycin I-type complex. 


This investigation was supported by a research grant from the United States 
Public Health Service and by funds supplied by the Rutgers Research and En- 
dowment Foundation. 


1 C.F. Dalgliesh, A. W. Johnson, A. R. Todd, and L. C, Vining, J. Chem. Soc., p. 2946, 1950. 

2H. Brockmann, Angew. Chem., 66, 1, 1954. 

* H. Brockmann and K. Vohwinkel, Naturwissenschaften, 41, 257, 1954. 

4S. A. Waksman, Antibiotics and Chemotherapy, 4, 502, 1954; 5, 409, 1954. 

‘H. Brockmann and N. Pfennig, Z. physiol. Chem., 292, 77, 1953; H. Brockmann, H. Linge, 
and H. Groéne, Naturwissenschaften, 40, 224, 1953; H. Brockmann, and H. Gréne, Chem. Ber., 87, 
1036, 1954. : 

®R. A. Manaker, F. J. Gregory, L. C. Vining, and 8. A. Waksman, Antibiotics Annual, 1954-65 
(New York: Medical Encyclopedia, Inc., 1955), p. 853. 

7L. C. Vining, personal communication. 

8 W. Braun, Bacterial Genetics (Philadelphia: W. B. Saunders Co., 1953). 

°S. A. Waksman and H. B. Woodruff, Proc. Soc. Exptl. Biol. Med., 45, 609, 1940; J. Bacteriol., 
42, 231, 1941, 


ON THE MECHANISM OF HYDROGENASE ACTION 
By Harry D. Peck, Jr.,* ANTHONY SAN Prerro,t AND Howarp Gestrt 


DEPARTMENT OF MICROBIOLOGY, SCHOOL OF MEDICINE, WESTERN RESERVE UNIVERSITY, CLEVELAND 
AND McCOLLUM-PRATT INSTITUTE, THE JOHNS HOPKINS UNIVERSITY, BALTIMORE 


Communicated by H. G. Wood, November 28, 1955 


The biological activation of molecular hydrogen to serve as an electron or hy- 
drogen donor is catalyzed by the enzyme hydrogenase. This enzyme, which 
occurs in a large variety of microérganisms, also appears to be the terminal com- 
ponent of the multienzyme systems responsible for photochemical formation of 
H, and for evolution of this gas in fermentations of carbohydrates, amino acids, 
purines, and related substances.! Both types of activity, i.e., oxidation and 
formation of hydrogen, have commonly been considered to be due to the same 
enzyme, particularly because certain organisms which contain hydrogenase 
behave as a reversible hydrogen electrode according to the following equation :? 


H, = 2H+ + 2e. (1) 


Recent investigations have demonstrated that electrons or hydrogen derived 
from activation of H» can readily be funneled into the well-known systems of 
electron transfer in which pyridine nucleotides, flavin, and cytochrome enzymes 
are prominent.* ‘ For purposes of purification and for study of the properties 
of hydrogenase itself, however, a dye such as methylene blue is ordinarily employed 


as the electron acceptor in the enzyme assay. The over-all reaction in this case is 
H, + methylene blue ~ methylene blue-2H. (2) 


An alternative assay system developed recently® is based on formation of molec- 
ular hydrogen from the reduced form of the ‘one-electron’? dye methyl viologen 
= —0.446 volts), Le., 


2 viologense + 2H* H, + 2 viologen. (3) 


A survey of hydrogenase preparations from a variety of organisms using the 
foregoing assays (reactions [2] and [3]) disclosed that both activities were usually 
present. Of particular significance, however, was the observation that certain 
hydrogenases appeared to catalyze only reaction (2).6 These results suggested 
that hydrogenase may be a ‘‘double-headed”’ enzyme with two prosthetic groups. 
According to this interpretation, one prosthetic group might be considered to 
mediate electron transport. to dyes or physiological acceptors of relatively high 
redox potential, while the other might be concerned with reduction of “one-electron” 
acceptors of low potential and also with formation of molecular hydrogen. 

In addition to the reactions already noted, certain hydrogenase preparations 
catalyze an isotope exchange reaction between water and molecular hydrogen as 
in equation (4) 


HDO + H. = + HD. (4) 
Closely related to this activity is the ability of such preparations to catalyze the 
ortho-para hydrogen conversion reaction.” 
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Reaction (4) has frequently been assumed to represent a “primary” activity 
of hydrogenase (i.e., the immediate and direct reaction catalyzed by the enzyme 
itself). The data available from several recent reports, however, suggest that the 
velocity of the exchange reaction in certain instances (in Azotobacter vinelandit® 
and Rhodospirillum rubrum’) is very slow as compared with “hydrogenation” 
activity (i.e., reaction [2]). This apparent discrepancy indicated the desirability 
of a quantitative comparison of hydrogenase activities in various preparations as 
measured by reactions (2), (3), and (4). Recent developments in the elucidation 
of metallo-flavoprotein oxidation-reduction enzymes” also suggested the need for 
further study of the metal components of hydrogenase. The results of such 
experiments and their implications for the mechanism of hydrogenase action are 
reported in this communication. 


EXPERIMENTAL 

The crude cell-free extracts were made as described by Peck and Cest,® except 
that the Escherichia coli extract was prepared by disruption of the cells in a 
Hughes press.!! Fraction 1 of Micrococcus lactilyticus was obtained by precipita- 
tion with ammonium sulfate (for details see Peck and Gest®); fraction 2 was 
derived from the same extract by further addition of ammonium sulfate to a 
saturation of approximately 40 per cent. 

Evolution of hydrogen from reduced methyl viologen (reaction [3]) and re- 
duction of methylene blue (reaction [2]) were measured essentially as described 
previously;> in the methylene blue assay, the center well of the Warburg vessel 
contained alkaline pyrogallol to remove residual traces of oxygen, and the enzyme- 
plus-buffer solution was supplemented with 10 mg. of crystalline albumin (Armour) 
per 1.2 ml. 

Exchange-activity assays (reaction [4]) were conducted in vessels of the design 
described by Hoberman and Rittenberg.!? The flasks were shaken at 25° C. at 
a rate of 130-140 strokes per minute on a reciprocal shaker with a 0.75-inch excur- 
sion. The liquid phase consisted of 10 per cent D,O; 0.05 M phosphate buffer, 
pH 6.5; enzyme preparation; and distilled water (final volume, 5 ml.). Except 
for control experiments, the actual quantity of each extract tested was that showing 
an activity of approximately 2,500-3,000 ul. He. utilized per hour with methylene 
blue as the electron acceptor. Just before each assay, 1.3 mg. of solid sodium 
hydrosulfite was added (in order to deoxygenate hydrogenase!) to each flask; 
the latter was immediately evacuated and flushed twice with tank hydrogen and 
then filled with hydrogen to a pressure of 0.5 atmosphere. The change in deu- 
terium content of the gas phase was measured at various intervals, using a mass 
spectrometer. 


RESULTS AND DISCUSSION 


All the preparations tested (with hydrosulfite present) showed linear exchange 
rates over a period of approximately 90 minutes. In the case of Clostridium 
butylicum the linear per cent D versus time curve did not pass through the origin, 
apparently because a small quantity of isotopic gas was released initially upon 
addition of hydrosulfite to the reaction mixture. In control experiments using 
the F. coli extract, it was established that the rate of exchange was proportional to 


— 
ast 


Vou, 42, 1956 BIOCHEMISTRY: PECK ET AL. 15 


enyzme concentration under the conditions employed and therefore was not limited 
by the rate of hydrogen diffusion into the liquid phase. 

a) Relationship of Exchange to Other Hydrogenase Activities——The exchange 
activities observed in the different microbial preparations and the corresponding 
activities determined with the other two methods of assay (reactions [2] and [3]) 
are listed in Table 1. 


TABLE 1 


Isorore Excuance, HyproGen Evoiution, AND Repuction 
AcTIvITIES IN DIFFERENT MICROBIAL PREPARATIONS 
Isotope Exchange 
Assay: Increase Methylene Blue 
in Atoms Per Cent Reduction Assay Evolution Assay 
Preparation Excess D/Hr/MI1 (ul He Utilized /Hr/Ml) (ul He Evolved/Hr/M1) 
M. lactilyticus: 
Crude 9.30 
Fraction 1 0.18 
Fraction 2 0.93 
. pasteurianum 10.60 
butylicum 1.00 
coli 1.80 
A. vinelandii 0 
. rubrum 0 


Thus methylene blue is readily reduced by certain preparations which do not 
show appreciable exchange activity, i.e., by M. lactilyticus (fraction 1), A. vine- 
landii, and R. rubrum.‘* It is significant that these particular hydrogenases also 
show essentially no activity by the evolution assay (reaction [3]). Rapid ex- 
change was observed only in preparations which showed activities in both the 


methylene blue reduction and the evolution assays. 

It may be noted that, in systems showing exchange, the ratio of evolution to 
reduction activity varies from approximately 2 to 4, except in the case of C. butyli- 
cum. The very large ratio observed with the latter preparation is due to the 
extreme sensitivity of this particular hydrogenase to inactivation by molecular 
oxygen (largely reversible) and other oxidizing agents. Thus, in contrast to 
other hydrogenases, reduction of methylene blue by C. butylicum extract can be 
successfully demonstrated only if very special precautions are taken to insure 
completely oxygen-free conditions. Molecular hydrogen is rapidly oxidized 
by this preparation with the one-electron dye benzyl viologen (29’ = —0.359) 
acting as the acceptor, but again only if a low redox potential is established In 
the evolution assay, the powerful reducing mixture of hydrosulfite plus methyl 
viologen apparently reverses inactivation, due to traces of oxygen, which may 
occur during experimental manipulations. From this standpoint, the evolution- 
assay procedure offers a distinct. experimental advantage, particularly since the 
sensitivity of most hydrogenases to inactivation by oxygen increases markedly 
with purification (see, e.g., Sadana and Jagannathan"). The conditions used in 
measuring the exchange reaction are not so strongly reducing as in the evolution 
assay, and this may account for the relatively low exchange activity shown by 
C. butylicum hydrogenase (i.e., the value listed is a minimum). 

The virtual absence of exchange activity in certain preparations capable of 
activating H. for reduction of dyes suggests modification of current views on the 
nature of bacterial hydrogenase and on the significance of the isotope exchange 
reaction as a primary assay for the enzyme (see below). 
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b) Metal Components of Hydrogenase.-—Recent experiments of Shug, Wilson, 
Green, and Mahler* have shown that MoO ; and flavin adenine dinucleotide are 
required for reduction of cytochrome ¢ by a purified hydrogenase from C. pas- 
teurianum. Similar results with M. lactilyticus hydrogenase have been obtained 
by Whiteley and Ordal,"* who studied H, oxidation with benzyl viologen acting 
as acceptor. 

There is compelling evidence indicating the presence of a second metal, viz., 
iron, as a functional component of hydrogenase. Of particular significance in 
this connection is the fact that hydrogenase action is markedly inhibited by CO. 
This is illustrated in Table 2, which represents a compilation of data obtained by 
a number of investigators in studies with a variety of hydrogenase preparations. 


TABLE 2 


CarBON MonoxipeE INHIBITION OF HYDROGENASE ACTIVITIES 
Volume 
Type of Activity Per Cent Per Cent Light 
Organism Preparation Measured co Inhibition Reversal Reference 

. pasteurianum Purified enzyme Methylene blue 20 75 None Shug and Wilson* 
reduction 

. butylicum Purified enzymef He evolution from 100 or Peck and Gestt 
methyl viologen 

. butyricum Intact cells He evolution from 
pyruvate 


Kempner and 
Observed Kubowitz;t 
Kubowitz§ 


None Joklikl| 


ooo 


. colt Cell-free extract Methylene blue 
reduction 
Intact cells He evolution from 


Stephenson and 


formate Stickland¢ 
Hoberman and 
Rittenberg** 
Hyndman, Burris 
reduction f and Wilsontt 


. vulgaris | Intact cells Exchange reaction Observed 


. vinelandii Cell-free extract Ferricyanide 
None 


* In Symposium on Inorganic Nitrogen Metabolism (Baltimore: Johns Hopkins University Press) (in press). 
t The specific activity of the purified enzyme was 2.5 X 108 ul He evolved/hr/mg protein nitrogen at 30° C. 
The purification procedure and properties of the enzyme will be described in detail elsewhere. 
Biochem. Z., 265, 245, 1933. 
Biochem. Z., 274, 285, 1934. 
Australian J. Exptl. Biol. Med. Sci., 28, 331, 1950. 
# Biochem. J., 26, 712, 1932. 
** J. Biol. Chem., 147, 211, 1943. 
tt J. Bacteriol., 65, 522, 1953. 


Light-reversibility of the CO inhibition, a characteristic of iron enzymes, has 
been reported in two instances, viz., for the production of Hz by intact cells of 
C. butyricum™ and for the exchange reaction in intact Proteus vulgaris." Failure 
to detect light-reversibility in a number of other cases does not appear to represent 
a critical argument against implication of iron, since in these instances dyes and 
other light-absorbing substances were generally present in high concentration and 
could well have interfered by absorption of effective wave lengths. 

Kempner and Kubowitz” obtained an action spectrum of the “hydrogen- 
evolving enzyme”’ in intact cells of C. butyricum by determining the efficiencies of 
CO inhibition reversal with light of different wave lengths. This spectrum has 
been interpreted by Warburg’ to indicate an iron enzyme, but not one of the 
porphyrin type. Since the activity of purified hydrogenase from C. butylicum 
(believed to be the same organism as (’. butyricum) is completely suppressed by 
CO, it appears that the inhibition studied by Kempner and Kubowitz was un- 
doubtedly due to combination of CO with the hydrogenase enzyme rather than 


* 
100 } 
{ 2s 
100 
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to inhibition of other enzymes concerned with the fermentation of pyruvate. 
Similarly, the sensitivity of /. coli hydrogenase to CO probably also accounts for 
the CO inhibition of Hs production from formate by the hydrogenlyase enzyme 
complex. 

Stimulatory effects of Fe+*+ on hydrogenase activity in partially resolved prepa- 
rations from W. lactilyticus® and C. butylicum (Peck and Gest, unpublished) have 
been noted. On the basis of all of these findings, together with those of early 
nutrition experiments! which indicated an iron requirement for the formation of 
the enzyme, it is concluded that iron is a functional component of hydrogenase. 

c) A Suggested Formulation of Hydrogenase Action.—On the basis of experiments 
with P. vulgaris, Krasna and Rittenberg’ have suggested that both the exchange 
reaction and the ortho-para interconversion are catalyzed by hydrogenase according 
to the following primary” reaction: 


H: + E=E:H- + Ht, (5) 


where E represents the hydrogenase enzyme. In the presence of D.O (i.e., D*), 
HD would be produced by reversal of reaction (5). It may be noted that the 
hydrogen metabolism of P. vulgaris closely resembles that of FE. coli, and, in con- 
formity with the results already described, extracts of the former organism also 
show both methylene blue reduction and evolution activities.® 

According to Krasna and Rittenberg,’ ‘any organism which possesses hydro- 
genase activity should be capable of carrying out both the conversion of para 
hydrogen and the exchange reaction.” The present results observed with the 
preparations from M. lactilyticus (fraction 1), A. vinelandii, and R. rubrum" are 
not consistent with their postulated mechanism. 

A consideration of the available data suggests the following formulation for 
the mechanism of action of hydrogenase: 


H, + E(Fe+*+) = E(Fe**) = E(Fe**)-flavin: 2H 
flavin Mo 
methylene blue 2 viologense + 2H* 


According to this scheme,”® molecular hydrogen is activated by a ferrous enzyme 
to form a reduced enzyme complex. The latter can be oxidized through the agency 
of a flavin component by methylene blue and similar acceptors. Alternatively, 
if molybdenum is properly combined with the enzyme, flavin: 2H can be reoxidized 
by one-electron acceptors such as viologen dyes. 

The reduction of viologens by Hy and formation of hydrogen from the reduced 
methyl dye would require the complete system, i.e., flavin and the two metals. 
Similarly, isotope exchange activity (i.e., conversion of D+ to HD) would be 
catalyzed only by the complete system. This formulation accounts for the corre- 
lation reported here between exchange and evolution activities. 

If the specifie molybdenum component were absent, the hydrogenase would be 
capable of reacting only with methylene blue, i.e., exchange and evolution activities 
would not be manifested. Preparations in this category would be those from 
R. rubrum, A. vinelandii, and M. lactilyticus (fraction 1). Further support for 
this view is found in the fact that molybdenum is required for reduction of cyto- 
chrome ¢ by resolved preparations of C. pasteurianum hydrogenase, whereas addi- 
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tion of this metal does not affect methylene blue reduction.‘ Similarly, addition 
of MoO; to partially resolved preparations of M. lactilyticus hydrogenase affects 
the rate of benzyl viologen reduction markedly but has relatively little influence 
on methylene blue reduction. 

Since all systems which produce hydrogen from physiological intermediates 
also show evolution activity, it would appear that a metal component permitting 
one-electron transfers is essential in normal H, formation by microérganisms. It 
is assumed that a low redox potential is required for this process (also for exchange 
activity), and it is possible that molybdenum may be of particular significance 
in this regard. Thus it has been suggested that the low /o’ value (—0.45 volt) 
of the iron-molybdoflavoprotein xanthine oxidase may reflect a contribution of 
molybdenum to the potential."2. The remarkable similarity of hydrogenase to 
milk xanthine oxidase has previously been noted.* ! 

The state of iron in hydrogenase is assumed to be ferrous for three reasons: 
(a) CO combines only with the ferrous form of iron enzymes;!> (b) cyanide, which 
specifically affects the ferric state, does not significantly inhibit hydrogenase 
activity under strictly anaerobic conditions (e.g., see Hyndman, Burris, and 
Wilson® and Hoberman and Rittenberg!*?); and (¢) the enzyme is very susceptible 
to inactivation by oxygen and other oxidizing agents.' *! Evidence recently 
furnished by Fisher ef al.'* indicates that the hydrogenase of P. vulgaris can be 
inactivated by oxygenation and that this type of inhibition is reversed by physical 
or chemical means. These considerations suggest that the initial activation of 
H, may occur by a reaction analogous to that described by Klotz and Klotz?* for 
the oxygenation of the iron protein hemerythrin. Further studies with purified 
hydrogenases will be required for a more detailed understanding of the initial 
reaction and the subsequent electron-transport pathways. 


SUMMARY 

The isotope exchange reaction between heavy water and molecular hydrogen 
has been measured in a number of different microbial hydrogenase preparations 
and has been compared with their ability (a) to reduce methylene blue with H, 
and (b) to produce H, from reduced methyl viologen. Certain preparations 
showing hydrogenase activity only by methylene blue reduction fail to catalyze 
the exchange reaction at an appreciable rate. Special importance is attributed 
to the correlation between absence of exchange activity and inability to evolve 
hydrogen from reduced methyl! viologen. Additional evidence implicating iron 
as an essential component of hydrogenase is cited, and a generalized formulation of 
the mechanism of action of the enzyme is presented. 
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PIGMENT CONVERSION IN THE FORMATIVE RESPONSES 
OF PLANTS TO RADIATION 
By 8. B. Henpricks, H. A. Borruwick, Anp R. J. Downs 
AGRICULTURAL RESEARCH SERVICE, PLANT INDUSTRY STATION, BELTSVILLE, MARYLAND 


Communicated November 25, 1955 


Many aspects of growth and maturation of plants are known from the close 
similarity of their action spectra to be controlled by a reversible photoreaction.' 
Among these responses are germination of seeds of some higher plants® and of fern 
spores,’® photoperiodic control of flowering,‘ coloration of some fruits,’ and etiola- 
tion.6 The photoreaction can be written as follows: 


Red radiation 
(Max. near 
6600 A) 
Pigment; + (Areactant) == Pigment, + (Changed reactant) 
(red-absorbing) (far-red 
absorbing) 


Far-red radiation 
(Max. near 7350 A) 


Darkness 
The essential contribution herein is the expression of the physiological response 
as a function of the fractional conversion of the pigment into the effective form, 
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whichever that may be, and the finding of min mum values for the products of the 
respective absorption coefficients and the quantum efficiencies for conversion of 
the two pigment forms in vivo. These ends are attained without knowledge of the 
chemical nature of the pigment or of its concentration. The method depends upon 
physiological information and upon the fact that the photoreaction is reversible 
and follows first-order kinetics with respect to energy in both directions. That 
the reaction rate is not limited by molecular collision in either direction is shown by 
temperature coefficients of unity for the photoreactions in both directions for 
germination of lettuce seed between 6° and 26° C.7 

The method of analysis is similar to the one used by Warburg and Negelein‘’ 
to find the products of the absorption coefficients and quantum efficiencies for 
dissociation of cytochrome oxidase-carbon monoxide by photoreversal in vivo of 
the poisoning with carbon monoxide of the respiration function of the oxidase. 

After irradiation with energy F in an absorption region A; + Adj, an unknown 
fraction F of the pigment in one form is converted to the other according to the 
differential equation dF/dE = k(1 — F). The solution of this equation is kE = 
log [1/(1 — F)]. If the energies required to give two physiological responses, | and 
2, have the relationship a = then log [1/(1 — = log [1/1 — F2)]. 
The same responses in the reverse direction, obtained with the energies of the re- 
versing radiation in an absorption region A, + Ads, will have a different ratio 8, 
with 8 log [1/F2] = log [1/F,]. Thus, if a and g are measured for any two arbi- 
trary degrees of physiological expression, values of fF; and F, can be calculated. 
The only limitation for the equations as set up is that the dark reaction be slow and 
the radiant energy used for one response be relatively free of that inducing the 
other, but AX need not be small. 

The products of the absorption coefficients, w (em?/gram molecular weight), 
and the quantum efficiencies, ¢, fc: both the red and the far-red form of the pigment 
can be determined from the fraction, F, of the pigment converted in a given spectral 
region by a given incident energy F (in einsteins/em?). Thus w6X = number of 
effective quanta = F'/E, where X is the fraction of the incident energy reaching the 
region of absorption. The molar extinction coefficient, in American usage, is 
[w/ 1,000] log, 10. 

Internode elongation of the bush Pinto bean is the first example considered. If 
these beans are grown with daily light periods of several hours under light sources 
such as “white” fluorescent lamps having low intensities in regions beyond 7000 A, 
the developing internodes are short. But if they are exposed after the end of each 
light period on four successive days to moderate energies, which can be delivered 
in a few minutes or less, in the region beyond 7000 A the internodes increase in 
length several fold; essentially the beans are changed from “bush” to “pole” at this 
stage of growth. If the plants are irradiated with red light following exposure to 
high energy in the far-red region of the spectrum at the end of the light period, the 
lengthening stimulus of the far red is nullified. The internodes may attain lengths, 
depending on the red energies, no greater than those of control plants. This is the 
reversal. Essentially the length of the second internode of the bean is a function of 
the fraction of the pigment converted by radiation, and it can be used as an index 
of the conversion irrespective of the nature of the function. 

Results of an experiment are given in Table 1. The mean internode lengths 
are the averages for lots of 24 plants. If two responses are taken as the internode 
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lengths 40.3 and 65.2 mm., respectively, the value of 8 for far-red radiation is 8.00/ 
4.00 = 2.00. From interpolation on a plot of the third columns against the first 
columns, the reverse times to give the two internode lengths with the arbitrary red 


TABLE 1 


Mean LENGTHS OF SECOND INTERNODES OF PINTO BEANS RESULTING FROM IRRADIATIONS FOR 
Various PERIODS WITH FAR-RED RAbDIATION AND WITH RED RapiaTION FOLLOWING SaTURATING 
Far-RED ExposuRES 


PERIODS OF MEAN INTERNODE LENGTHS PERIODS OF MEAN INTERNODE LENGTHS 
IRRADIATION FoLLOWING IRRADIATION (Mm.) IRRADIATION Fo IRRADIATION (M™m.) 
(Mrn.)* Far-red Red (Mrn.)* F d Red 
2.00 
4.00 
8.00 
16.00 


* Each radiation source had a power equivalent to about 0.25 milliwatt/em? at the 
wave length of maximum action integrated over the region of action. 


radiation source are 5.95 and 1.63 minutes, respectively, giving a = 5.95/1.63 = 
3.65. The corresponding solutions of F; and F», are 0.10 and 0.316. A measure 
of precision is indicated by calculation for 8 = 2.00 and a = 3.30 or 4.00. The 
respective values of F; and F, are 0.12 and 0.348 and 0.08 and 0.284. The value of 
k for the far-red irradiance is accordingly 0.35 per minute of irradiation. The com- 
plete curve of physiological response as a function of pigment conversion by far- 
red radiation as calculated for this value of k is shown in Figure 1. The fact that 


100 
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Length of internode in mm. 


ie) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Fraction of pigment in far-red absorbing form 


Fic. 1.—Internode lengths of Pinto beans as a function of pigment conver- 
sion: circles are for irradiation in the red part of the spectrum of beans with 
the pigment initially in the red-absorbing form, and crosses are for irradiation 
in the far red with the pigment initially in the far-red-absorbing form. 


the points for both red and far-red radiation fall along the same curve validates the 


first-order kinetics. 
Values of w6X are next determined for the pigments in the spectral regions of 
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maximum action. These values are based upon internode lengths for beans from 
which one of the two primary leaves was removed so that the plants might be ir- 
radiated in a spectrum on four successive nights following their growth periods 
under fluorescent lights. One group of 18 plants was irradiated in the region 7100— 
7600 A and another group in the region 6300-6700 A after they had received a satu- 
rating exposure to far-red radiation. The first group received an incident energy of 
9.1 X 10-7 einstein/em? in each daily irradiation and the second 1.84 & 1077 ein- 
stein,em?. The fractions of the pigment converted for the two groups were 0.81 
and 0.21, respectively, as based on the variation of internode lengths of similar 
beans, with a single leaf, as a function of incident energy. The values of wX are 
0.09 107 and 0.11 X 10° for the far-red and the red-absorbing forms of the pig- 
ment, respectively. 

Before this procedure is extended to other examples, attention is again directed 
to the fact that the “mean”’ physiological response of bean internodes serves as an 
index of the pigment conversion. Each internode value listed in Table 1, for ex- 
ample, is a mean obtained from 24 individuals that varied in internode lengths. 
Reasons for the variations were many, including variation in the fraction of the 
pigment converted by a given incident energy. 

The procedure used for analysis of internode growth may be applied with equal 
validity to germination of seeds, despite differences in the character of results ob- 
tained in the two cases. The length of an internode following a light treatment 
may have any value between that for the untreated and that for saturated control 
plants, but the response of a seed following light treatment is either germination 
or lack of germination. After red irradiation, some imbibed seeds germinate upon 
being held in darkness at a favorable temperature while others remain dormant. 
At a slightly higher red energy the germination requirements of these same seeds 
and those of an increment are met. The seeds of this increment are different from 
all others in that their energy requirements for germination are greater than the 
Srst energy given but less than the second. The crucial point is that the analysis 
is based only on such an increment, and this increment is the same whether reached 
by promotion of germination with red or inhibition of germination with far-red 
radiant energy. 

The germination responses of lettuce, Lactuca sativa,’ and of pepper grass, 
Lepidium virginicum,® seeds to an increasing series of red and far-red energies have 
been measured, as illustrated for L. virginicum in Figure 2. Action spectra giving 
the incident energies required for given percentage germinations are also known.” ® 
All the information thus is available for calculations of the germination response as 
a function of the fractional conversion of the pigment and for finding the product 
of the absorption coefficient, the quantum efficiency, and X, that is, eX. The 
germination functions for the two kinds of seed as calculated are shown in Figure 3. 

Lepidium seed imbibed in 0.2 per cent KNO; solution requires about 4.5 & 105 
ergs/cm? at the action maximum of 6400-6600 A for 50 per cent germination at a 
constant temperature of 20°C. This per cent germination requires 0.17 conversion 
of the pigment. An energy of 7.0 X 10‘ ergs/cm? at the action maximum of 7200- 
7500 A is required for 50 per cent inhibition of fully promoted seed. The values 
of weX calculated from these quantities are 0.067 107 for the red-absorbing form 
of the pigment and 1.9 X 10° for the far-red-absorbing form. For lettuce seeds 
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these values are 5.0 X 10’ and 0.11 X 107, respectively, with 50 per cent germi- 


nation corresponding to 50 per cent pigment conversion. 
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Fic. 2.—Variation in germination of Lepidium virginicum seeds with times 
of irradiation for promotion in the red part of the spectrum and with inhibition 
in the far-red after full promotion. Conditions under which the seeds were 


held are indicated. 
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Fic. 3.—Variation in germination of Lepidium virginicum seeds and of Grand 
Rapids lettuce seeds with pigment conversion. The circles represent the 


interpolated values for 50 per cent germination. 


An immediately evident fact is that the values of w/X vary for the several ob- 
jects, but the question whether the variation is in w, ¢, or X involves further con- 
siderations. These considerations can be based on results for germination of L. 
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virginicum seed as reproduced in Figure 2. As conditions affecting germination are 
varied, the several lines expressing promotion of germination as a function of inci- 
dent energy are displaced without change of slope. This means that the fraction of 
pigment converted for a given germination is constant but that the incident energy 
required for this conversion varies. Since the seed is not changed in its light- 
absorbing or light-scattering characteristics by the several conditions, XY is approxi- 
mately constant. The conclusion is that ¢, the quantum efficiency for pigment 
conversion, varies with conditions of seed germination. A deduction to account 
for variation in ¢ is that a reactant other than the pigment is involved in the photo- 
reaction in both directions. 

A probable maximum value of the quantum efficiency for pigment conversion 
per chromophoric group is 1.0. The value of X also cannot be greater than 1.0, 
and an order of magnitude for seeds and leaves would be in the range of 0.1-0.5 or 
less. The maximum values found for w@X are 5.0 * 107 for red absorption in 
lettuce seed and 1.9 X 107 for far-red absorption in Lepidium seed. The minimum 
value of the molecular absorption coefficients of the pigment, w/1,000, in both 
forms at the action maxima, then, are of the order of 10°. The oscillator strengths, 
accordingly, are of the order of 1.0, the intrinsic time constant of the excited states 
being of the order of 10~° second. 

The values of absorption coefficients of the two pigments in the region between 
4000 and 5200 A are of the order of 0.01 times those at the action maxima in the 
red and far-red portions of the spectrum. This was shown’ by the rate of approach 
with increasing radiant energy to an equilibrium physiological response in the re- 
gion 4000-5200 A of objects in which the pigment systems were displaced to the 
two extremes by absorption in the red or far red. The method is independent of 
possible overlapping absorptions. The significance here is that cyclic tetrapyrolle 
systems, such as are present in porphyrins and might naturally be suspected here, 
are eliminated as possibilities for the pigment structures on account of their strong 
absorptions in the region in question. Instead, the pigments are probably open- 
chain conjugate systems, which include the phycocyanins of the blue-green algae 
as possibilities. The molecular extinction coefficient of c-phycocyanin at the ab- 
sorption maximum of 6150 A is 2.76 & 10° as measured by Svedberg and Katsurai.'° 
The corresponding probable value of the molecular absorption coefficient referred 
to a single chromophoric group is 1.59 X 10°, based on four groups in each molecule. 
Allophycocyanin has an absorption maximum at 6500 A!! and its relative absorp- 
tion coefficients at other wave lengths correspond closely to those of the red-ab- 
sorbing form of the pigment. Polyenes, of types as yet unknown, are not excluded 
as possible structures of the pigments. 

The details of the photochemical reaction are of interest. With the reactant’s 
presence now demonstrated, its general nature can be anticipated. The photo- 
chemical reaction probably involves the radiationless transition of the molecule 
from an excited singlet- to a triplet-state. The triplet-state free radical probably 
accepts hydrogen from the reactant in one form and donates hydrogen to the 
oxidized reactant in the other form. The variations in quantum efficiency arise 
from the variations in the ratio of the oxidized and reduced forms of the reactant 
in the biological object depending upon ambient conditions. For the photoreaction 
to be first order with respect to energy in the presence of the reactant implies con- 
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tinuous association of either the oxidized or the reduced reactant and the pigment, 
as is characteristic of condensed systems. 

The pigment in its far-red-absorbing form is evidently biologically effective. 
This deduction rests on three arguments, none of which is rigorous. The first 
argument is that, in two of the three objects considered, half-maximum physiologi- 
cal response was obtained with conversion of the order of 0.1 of the pigment from 
the red- to the far-red-absorbing form, and conversion of 0.75 approached physi- 
ological saturation. This type of dependence upon the fraction of the active far- 
red-absorbing form corresponds to that usually found in enzymatic control of a 
reaction. The second argument rests upon the fact that some seeds with the pig- 
ment in the red-absorbing form will lie dormant for about a century but will 
germinate in a few days when exposed to light, with consequent appearance of the 
far-red-absorbing pigment form. Inhibitors of germination present in seeds 
usually disappear after a year or more, so that discharge of the final dormai.cy 
would seem to depend upon the far-red-absorbing pigment form acting as an acti- 
vator rather than the red-absorbing form acting as a suppressor of an inhibitor. 
The third argument is that photoreversibility in photoperiodism is lost, an indica- 
tion of biological change, shortly after the pigment is converted to the far-red- 
absorbing form. 

The pigment is not sensitive to destruction by photo-oxidation in either of the 
two forms. It is freely available for photoconversion. This follows from the fact 
that the pigment system was fully effective in the green bean leaf after an 8-hour 
light period at high energy, as shown by the results given in Table 1. During the 
day the fraction of the pigment in the far-red-absorbing form in plants growing in 
the open is determined by the ratio of the effective red and far-red radiation in the 
radiation from the sun and by the ratio of the oxidized to the reduced form of the 
reactant in the leaf. With darkness the far-red form of the pigment changes 
slowly by a thermal reaction toward an equilibrium with the red-absorbing 
form.” 

Summary.—Variations of physiological response with the fraction of pigment 
conversion as obtained from the reversibility and the first-order character of the 
conversion with respect to energy in the controlling reversible photoreaction were 
measured for internode lengthening and seed germination. Knowledge of this 
fraction permitted calculation of absorption coefficients and quantum efficiencies 
for pigment conversion. The method is a general one applicable to photoperiodism 
and other responses controiled by the reversible photoreaction. 
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PROGRESSIVE CONFIGURATIONAL CHANGES, DURING CELL 
DIVISION, OF CELLS WITHIN THE APICAL MERISTEM* 


By Epwin B. Matzke 
DEPARTMENT OF BOTANY, COLUMBIA UNIVERSITY 
Communicated by Franz Schrader, October 20, 1955 


During the past several decades the three-dimensional shapes of cells which are 
essentially undifferentiated morphologically have been studied by a number of in- 
vestigators. Many different organisms have been used. Some research has also 
been done on configurational changes, during cell division, in epidermal layers. 
Thus far, however, no analysis has been made anywhere of the shape changes, step 
by step, which take place during the course of cell division. Furthermore, no study 
whatever, from this standpoint, has heretofore been carried out on cells deep within 
a tissue. In view of the increased emphasis, both morphological and physiological, 
placed upon terminal and lateral plant meristems, the internal cells, within the 
stem apex, were chosen for this study of the progressive shape changes of dividing 
cells, 

Previously comparisons have been drawn between the shapes of interphase cells 
and those of cells in division in the epidermis of the apical meristem of Anacharis 
densa (Matzke'); Mozingo? has traced the configurational alterations, during 
division, of epidermal cells of living roots of Phleum pratense. Various aspects of 
cell division have also been considered in a series of publications by Lewis.* 

It is now known as a result of the research of numerous workers that cells which 
are essentially undifferentiated morphologically, in tissues, have an average of ap- 
proximately 14 faces and that no one type is achieved. Much of the literature deal- 
ing with this has already been summarized (Matzke‘). Cells with 8 hexagonal 
and 6 quadrilateral faces, as in the tetrakaidecahedra of Kelvin,’ occasionally occur, 
but they are infrequent. In analyzing the factors responsible for shape determina- 
tion, experiments have been performed on compression of lead shot by Marvin’ and 
on bubbles in foam by Matzke.?- Comparisons have been made between the shapes 
of cells in tissues and grains in a metal alloy by Williams and Smith*® and by Smith.® 
Detailed analyses of similarities and differences in the data presented below and 
the results of previous investigators will be published elsewhere. 

Materials and Methods.—Vigorously growing vegetative stem apices of A. densa 
were fixed in Craf, sectioned lengthwise at 40 u, and stained in Feulgen and in 
ruthenium red. Only the median section of a tip was used, and it was examined 
under oil immersion with a Leitz Ortholux microscope. The cells were selected 
from the region of the apex above the visible origin of lateral primordia, and they 
were at least 4 cells from the periphery in any direction. This region is shown by 
the brace in Figure 1. 
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The 400 interphase cells tabulated below were taken from 193 stem tips, and the 
400 dividing cells came from 329 apices. Some of the sections were utilized for 
both the dividing and nondividing cells. Since the 
slides could not be tilted under the microscope, the 
configuration of the individual cells was determined 
by careful and prolonged focusing up and down. 
Only those cells that could be unmistakably worked 

out in three dimensions were utilized, and con- 
sequently each median section yielded only one or a 
few cells that were clear and could be recorded. A 
freehand sketch was made of each cell used, showing 
all the contacts. Some of the cells were also illus- 
trated with camera lucida drawings, eight of which 
\ are shown in Figures 2-9. These drawings are 
\y ZB slightly schematized, and minor curvatures of the 

\ walls are not reproduced. 
_ Fic, 1.—Camera lucida draw- The dividing cells were placed in seven groups 
ing of a longitudinal section of a 

stem tip of A. densa. The cells according to the stage of division (Table 1). In 
studied are in the region indicated early prophase (Fig. 3) the Feulgen stain was more 

by the brace. X68. 
’ intense than in interphase, and the chromatic 
material had become more evident. By middle prophase (Fig. 4) the chromosomes 
were more readily distinguishable, while in late prophase (prometaphase) the 
chromosomes were shorter and more deeply stained and the nucleoli and nuclear 


TABLE 1 


NuMBER OF Faces 1n 400 INTERPHASE CrELLs, IN 400 Divipine CELLS, IN VARIOUS 
SraGes or Division, AND IN 150 DAUGHTER CELLS, WITHIN THE APICAL MERISTEM 
(SHoot Apex) or Anacharis densa (ELopEA) 


Interphase 400 1 7 22 61 78 101 
Stage of mi- 

tosis: 
Prophase 
Early 
Middle 
Late 
Metaphase 
Anaphase 
Telophase 
Early 
Late 


Total di- 
viding 
Daughter 150 


* Standard error of the mean. 
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12.61 + 0.10 


membranes were no longer evident (Fig. 5). Metaphase and anaphase are illus- 
trated in Figures 6 and 7 and were determined in accordance with the usually 
accepted criteria. Early telophase cells (Fig. 8) varied from those in which the 
cell plate was barely visible to those in which it nearly but not quite came in con- 
tact with the original walls of the parent cell. In late telophase (Fig. 9) it made 
contact with those walls. 
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Data.—It is evident from Table 1 that the 400 interphase cells had an average 
of 13.85 faces and that the range extended from 9 to 21. Cells with 14, 13, 15, and 
12 contacts occurred most commonly, in that order of frequency. An interphase 
cell with 14 faces is illustrated in Figure 2. The interphase cells recorded were not 


3 


INTERPHASE EARLY PROPHASE 


4-4-6 (14) 6-3-4-3 (16) 


MIDDLE PROPHASE LATE PROPHASE 


3-7-5-1 (16) 5-5-4-3 (17) 


Figs. 2-5.—Camera lucida drawings, slightly schematized, of four characteristic cells, one in 
interphase and one each in early, middle, and late prophase. 2,500. The formula for the 
combination of faces is given below each cell, i.e., 4-4-6 is 4 quadrilaterals, 4 pentagons, 6 hexagons, 
and the total number of faces for each cell is in parentheses. The cells, from different, tips, are 
drawn in the same relative position, and the apex of the tips is indicated by the arrow. 


in contact with dividing cells. The kinds of faces varied from triangular to non- 
agonal (Table 2), with hexagons, pentagons, and quadrilaterals by far the most 
abundant. The average number of faces for 1,150 cells in diverse tissues, previ- 
ously studied by various investigators, was 13.86. These cells were presumably 
all in interphase. 
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METAPHASE ANAPHASE 


4-6-6-0-1 (17) 3-7-6-1 (17) 


EARLY TELOPHASE LATE TELOPHASE 


5-3-4-1 (13) 


Fics. 6-9.—Camera lucida drawings, slightly schematized, of four characteristic cells, one each 
in metaphase, anaphase, early telophase, and late telophase. 2,500. The formula for the 
combination of faces is given below each cell, and the total number of faces for each cell is in 
parentheses. In late telophase (Fig. 9) the formulas for the combinations of the two daughter 
cells into which the parent cell divided are also given. The cells, from different tips, are drawn 
in the same relative position, and the apex of the tips is indicated by the arrow, 
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The 400 cells referred to above occurred in 96 different combinations of faces, 
12 of which are listed in Table 3. The most common combination, represented 
by 38 cells, had 4 quadrilateral, 4 pentagonal, and 6 hexagonal faces (4-4-6)—-14 
faces in all. One of these cells is illustrated in Figure 2. Other common combina- 
tions were 4-4-5, which was found in 28 cells; 4-4-4, which occurred in 25 cells; 
and 3-6-4, which was represented by 22 cells. These 4 combinations were also 
the most common in the 1,150 undifferentiated cells previously studied by various 


TABLE 2 


Kinps oF Faces 400 INTERPHASE CELLS, IN 400 Drvipina CELLS, IN VARIOUS 
Sraces oF Division, AND IN 150 DauGHEeR CELLS, WITHIN THE 
APICAL MERISTEM OF Anacharis densa 
No. Tri- Quadri- Pen- Hex- Hep- Octag- Nonag- Dee- 
CeLLs angular lateral tagonal agonal tagonal onal onal agonal 


Interphase 400 5 1,605 — 1,837 1,844 239 10 1 0 


Stage of mitosis: 
Prophase 
Early 249 322 326 73 
Middle 306 418 394 101 
Late 99 152 120 38 
Metaphase K 228 356 338 67 
Anaphase 195 295 267 70 
Telophase 
Early 5 233 300 327 80 
Late 308 424 421 119 
Total dividing 400 1,618 2,267 2,193 548 
Daughter 150 597 691 506 91 


TABLE 3 


Most FREQUENTLY OCCURRING COMBINATIONS OF FACES IN 
INTERPHASE, Divipina, AND DAUGHTER CELLS WITHIN 
THE APICAL MERISTEM OF Anacharis densa 
COMBINATION ———No. ~ 


———No. Sipes PER Face———— 400 400 150 
6 7 Interphase Dividing Daughter 


1010 
128 
253 
14° 
224 
28? 
38! 
101 


ww 


* Superior numbers indicate order of frequency. 


authors. The 96 combinations of faces mentioned above were in 202 different 
arrangements, as determined by the construction of conventionalized models of 
each combination that was found more than once. In 27 arrangements there 
were 2 or more cells that were mirror images of each other Four of the 400 inter- 
phase cells had 6 quadrilateral and 8 hexagonal faces, as in the tetrakaidecahedra of 
Lord Kelvin. However, 20 other combinations occurred more often than the 
6-0-8, and 4 others, in addition to the 6-0-8, were found four times, making 25 in 
all that were represented by 4 or more cells. It is evident, therefore, that there 
was no one predominant “type” in these interphase cells; that they varied in num- 


Total 
15 0 0 986 
22 3 0 1,252 
6 1 0 417 
0 0 1,009 
9 0 837 
13 0 0 953 
7 2 1 1, 283 
89 6 1 6,737 
3 0 O 1,892 
Tora. 
Faces 
11 x 113 
12 16! 
12 94 
13 sh 6° 
14 7 44 
15 1 15! 
16 1 4 104 
16 1 2 11° 
16 2 3 105 cs 
17 1 3 113 


VoL. 42, 1956 BOTANY: E. B. MATZKE 31 


ber, kinds, combinations, and arrangements of the faces, but within the general 
limits outlined above; an‘ that only | per cent of the total had 6 quadrilateral and 
8 hexagonal faces, as in the Kelvin tetrakaidecahedra. 

The average of the 400 dividing cells was 16.84 faces (Table 1), the range extend- 
ing from 13 to 22. Cells with 17, 16, 18, and 15 contacts occurred, respectively, 
most frequently. Pentagons, hexagons, and quadrilaterals were the most numerous, 
but other kinds of faces—heptagonal, octagonal, nonagonal, and triangular—were 
present far more often in the dividing than in the interphase cells (Table 2). 

In Table 1 the number of faces is recorded for the stages of division. By early 
prophase the average had risen from the 13.85 of interphase to 16.43; in metaphase 
it was 16.82, and by telophase it was above 17. Thus the increase between inter- 
phase and early prophase was more than 2.5, while that between early prophase 
and telophase was definitely less than 1. 

The 400 dividing cells occurred in 170 combinations of faces, of which 96 were 
found only once. The 5 found most frequently are listed in Table 3. The com- 
bination 3-7-4-1 was found fifteen times, 4-5-6-1 eleven times, and 4-5-7-1 eleven 
times, while 3-7-5-1 (Fig. 4) and 4-6-4-2 each occurred ten times. Unlike the com- 
monest combinations of interphase and daughter cells, those in the dividing cells 
all had at least one heptagonal face. In these 170 combinations there were 366 
different arrangements of the faces, of which 343 were each represented by a single 
cell. There was one cell, an anaphase, which had 6 quadrilateral and 8 hexagonal 
faces, as in the tetrakaidecahedra of Kelvin. 

The 75 late telophase cells had cell plates which extended to the walls of the 
parent cells, so that the configurations of the 150 daughter cells thus formed could 
be determined. These 150 daughter cells had an average of 12.61 contacts, the 
range extending from 10 to 16 (Table 1). Twelve- and thirteen-faceted cells were 
most numerous. Pentagonal, quadrilateral, and hexagonal faces were more com- 
mon than any other kinds (Table 2). 

The 150 daughter cells showed 49 different combinations of faces, 24 of which 
were found just once, and 85 different arrangements of the faces. The combina- 
tions occurring most frequently (Table 3) were 4-4-4, 16 cells; 3-6-3, 14 cells; 
4-4-3, 11 cells; 5-2-5, 9 cells; 3-6-4, 6 cells; and 4-4-5, 6 cells. All of these had 
only quadrilateral, pentagonal, and hexagonal faces. These were also common 
combinations in the interphase but not in the dividing cells. 

The same daughter-cell combination was often derived from different parent- 
cel] combinations. Thus 13 different parent-cell combinations, varying in number 
of faces from 15 to 18, gave rise to the 16 daughter cells with the combination 4-4-4, 
and 14 different parent-cell combinations gave rise to the 14 daughter cells with the 
combination 3-6-3. The 2 daughter cells derived from the same parent cell usually 
had different combinations of faces and nearly but not quite the same number of 
faces. Two such daughter cells, one with 12 and the other with 13 faces, are il- 
lustrated in Figure 9. One of the 150 daughter cells tabulated had 6 quadrilateral 
and 8 hexagonal faces, as in the Kelvin figures. 

If the cell plate is hexagonal in outline, then the sum of the faces of the 2 daughter 
cells will exceed by 8 the number of faces of the parent cell, as in Figure 9; if the 
cell plate is pentagonal, the sum of the faces of the daughter cells will exceed those 
of the parent cell by 7; if heptagonal, the sum will be greater by 9, etc. Inspection 
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of Figure 9 shows further that in the division of a parent cell to form 2 daughter 
cells by means of a hexagonal cell plate, 6 faces are added to the contiguous cells, 
in addition to the 8 by which the sum of the faces of the 2 daughter cells exceeds the 
number of faces of the parent cell; in other words, division in this way results in 
the addition of 1 cell and a total of 14 faces to the cell mass. If the cell plate is 
pentagonal in outline, a total of 1 cell and 12 faces are added to the cell mass; if 
the cell plate is heptagonal, 16 faces are added to the cell mass. This occurs in- 
dependently of the number of faces of the parent cell. In the 75 late telophase 
cells studied, 48 had hexagonal cell plates, 18 had heptagonal cell plates, and in 9 
the cell plates were pentagonal. Therefore, a total of 14 faces was added to the 
cell mass of the growing apex as a result of the great majority of the cell divisions. 
This is in spite of the fact that 117 of the 150 daughter cells had fewer than 14 
faces (Table 1), and it is independent of the total number of faces of the late telo- 
phase cells, which ranged from 15 through 21. 

After interphase, dividing, and daughter cells (Table 1) had been recorded, a 
sampling was made to determine the ratio of dividing to interphase cells in the 
apical meristem. In the same limited region of the tip utilized for assembling the 
cells in Table 1, there were, in the median sections of 102 apices, 4,000 interphase 
cells and 322 cells in the various stages of division; 56 of these were in late telo- 
phase, and there were therefore 112 daughter cells. Utilizing the ratio of inter- 
phase, dividing, and daughter cells thus determined—including daughter rather 
than late telophase cells—and the averages established in Table 1, a weighted aver- 
age was computed for all the internal cells of the stem tip. This weighted average 
was 14.0. While the exact number 14.0 may not be important, there can be 
little doubt that the internal cells of the stem apex of Anacharis have an average of 
14, or of close to 14, faces when all are considered together. 

Conclusions and Summary.—From the foregoing data it may be concluded that 
interphase cells within the apical meristem have an average of slightly less than 14 
faces (13.85). No one type is achieved, since the 400 cells occurred in 96 different 
combinations of faces, some of which, however, were much commoner than others. 
There were 4 cells with 6 quadrilateral and 8 hexagonal faces, as in the tetrakaide- 
cahedra of Kelvin, but 38, for example, with 4 quadrilaterals, 4 pentagons, and 6 
hexagons (Table 3). Four hundred cells in various stages of division had an aver- 
age of 16.84 faces. Most of the increase in number of faces occurs either prior to or 
by early prophase; a slight further increase takes place during the stages of di- 
vision, and some additional contacts are also gained after the daughter cells have 
been formed. Swift! and several other investigators have reported that increase 
in DNA occurs in interphase preceding mitosis. The major increase in number of 
contacts in dividing cells is in essential agreement with this. 

There is no one combination of faces which predominates during any of the stages 
of division, the 400 dividing cells having oecurred in 170 combinations of faces. 

The average number of faces for 150 daughter cells was 12.61—fairly close to 
that of the interphase cells. The combinations which occurred frequently in the 
daughter cells were also commonly found in the interphase cells but not in the cells 
in division. 

Additional faces originate by division of contiguous cells; no evidence for “ gliding 
growth” was obtained in the apical meristem. 
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It may further be concluded that cell division results in a plexus of diversified 
three-dimensional forms which arise from the interaction of numerous factors but 
that this occurs within fairly rigid mathematical limits. Inorganic systems which 
have previously been studied—bubbles in foams, compressed lead shot, grains in 
a metal alloy—show nothing which is comparable with the progressive configura- 
tional changes which characterize the stages of cell division. This would be ex- 
pected, since cells originate in this intricate way from pre-existing cells. There is 
yo one configurational “type’’ which is achieved, either in interphase, in division, 
or in the daughter cells resulting from division. When the numerous factors which 
play a part are considered—surface forces, contact and pressure, elasticity of the 
cell walls, directions of cell division, volume relationships of the cells, unequal 
growth of the cells, necessities of spherical geometry, and inherited form tendencies 
—no one “type” would be expected. In spite of this, there are both regularity and 
order, as indicated above, and the variation and spread which do occur are within 
definite limits, as the results show. These limits are in general more circumscribed 
than those in inorganic systems. 


* This investigation was assisted by grants from the Eli Lilly Company and from the Fund for 
Research and the Higgins Fund of Columbia University. 
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THE OPTICAL ROTATORY POWER OF HELICAL MOLECULES* 
By Donaxp D. Firrst anp Joun G. Kirkwoop 
STERLING CHEMISTRY LABORATORY, YALE UNIVERSITY, NEW HAVEN, CONNECTICUT 
Communicated November 15, 1955 


Interest in the physical properties of helical molecules has increased in recent 
years, largely because of the helical configuration of many polypeptide chains. 
Since a helix has no elements of reflection symmetry, a molecule with a helical 
configuration is optically active even in the absence of asymmetric carbon atoms. 
In this article we shall apply the polarizability theory of optical activity to the 
determination of the specific rotation due to the helical configuration for an in- 
finitely long helical molecule. 

According to the polarizability theory of Kirkwood,! the specific rotation of 
plane-polarized light of sodium D wave length by an optically active compound is 
given by 
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.(n? + 
= 4.930 x l 
la]p x 3M (1) 
Here n is the index of refraction of the solution, and M is the molecular weight of 
the compound. For molecules having substituent groups with polarizability 
tensors of cylindrical symmetry, the molecular rotatory parameter g is given by 
N 
=1 
Gx — 3 
| Ri? 
Ri, : R.. 


The molecule is divided into groups each of which has a polarizability tensor of 
cylindrical symmetry, and the summation in equation (2) is taken over each pair 
of groups. R, is the vector drawn from the center of mass of the molecule to the 
center of mass of the 7th group, e is the unit vector along the axis of optical sym- 
metry of the group, @ is the mean polarizability of the group in cubic angstroms, 
and 8 is the anisotropy ratio of the group. 

The centers of mass of the atoms forming the main structural chain are placed 
on a right-handed cylindrical helix whose equations are x = a cos 6, y = a sin @, 
and z = b@. In equation (2) R is taken as the distance from the origin to a point 
on the helix. We shall consider here only cases where the polarizability tensors 
are symmetrical about the tangent to the helix. Therefore, the unit vector e 
becomes the unit tangent to the helix and is given by dR/ds, where ds is an element 
of are length. A case where the symmetry axis of the polarizability tensor is 
perpendicular to the helix has been discussed elsewhere.” 

In applying equations (1) and (2) to helical molecules, we assume that the 
polarizability is continuously distributed over the length of the helix. The sum- 
mation may then be replaced by an integration over the variable 6, so that 


OL OL 
“12 ao. | X ex) (3) 
0 0 


where a is the mean polarizability per unit radian and @ is the anisotropy ratio. 
Both a and 8 are assumed to be constant along the helix. The limits of integration 
are zero and 6,, the total angle of revolution. The factor '/2 is introduced because 
the double integral counts each pair of angle elements twice. 

The trigonometric functions in g® combine, so that the expression becomes a 
function of — 6, only. Letting y equal b/a, we find that 


OL OL 

g = an, f f (Ox) d6;, (4) 
9 sin @ — 2(1 — cos 6) 
[2(1 — cos 6) + 


= 


3(sin + 0)? | 
[cos 2(1 — cos + 


The double integral in equation (4) is easily replaced by 


| 
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0 Ok — OL 


which can be integrated by parts to yield 


OL 
— Ox) Gx) dix 

if f(0,) iseven. In the following discussion, the subscripts on the variable 6, will 

be dropped for simplicity. 

For an infinitely long helix, @; — @ can be approximated by just 6,, which equals 
2rv, where v is the number of turns of the helix. The mean polarizability per 
radian, a, can be replaced by a:/2x, where a is now the polarizability per turn. 
Unfortunately, f(@) cannot be integrated analytically. However, for large @, 
(9) may be approximated by 

[@ sin @ — 2(1 — cos 6) ](cos 6 — 2?) 


( 


which can be integrated analytically. The integral is therefore divided into two 
parts: 


= f(8) a+ «| (6) 


The first integral is evaluated by numerical integration. The second is integrated 
by parts until it reduces to the integral cosine, 


Ci(e) = dx, 


the value of which can be obtained from tables. 

Equation (6) is the most general expression for g“” for an infinitely long molecule 
whose polarizability tensors are symmetric about a uniform cylindrical helix. We 
shall now apply equation (6) to a typical alpha helix with 3.69 amino acid residues 
per turn. The average radius a is 1.81 A, and b is 5.44/27 A’ If ¢ is taken as 
67, the first integral has a value of 0.467 and the second Ci(127) — Ci(6x)/y' 
If € is greater than 6, the result is substantially the same. Therefore, for a right- 
handed alpha helix, 


g® = 1.15 XK 


56 2R2 2 2 
(7) 


3m 


where m is M/v, the molecular weight per turn. 

To estimate the magnitude of the contribution to the optical activity by a helical 
structure, we shall compute the theoretical rotation of a hypothetical alpha helix 
of polyglycine, the only polypeptide having no asymmetric carbon atoms. 

The mean polarizability for the various substituent groups of polyglycine is 
determined from the atomic refractions‘ by the relationship 
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(8) 


The values of a; in cubic angstroms are 1.81 for CH», 1.82 for C = O, and 1.42 
for NH. If we assume that the mean polarizabilities are additive, the mean 
polarizability per turn, a1, is 3.69 X 5.05 A. The anisotropy ratio 8 may be taken 
as about '/; for the glycine residue, and the molecular weight per turn is 3.69 x 
57.05. An estimated specific rotation for an infinitely long right-handed alpha 
helix of polyglycine in an aqueous solution with index of refraction 1.35 therefore 
is +132°. 

The available experimental measurements of the optical activity of polypeptides 
cannot be applied to test the validity of equation (7) because these polypeptides 
are mixtures of right- and left-handed helices. Some polypeptides chains con- 
taining constituents with asymmetric centers may, because of asymmetric syn- 
thesis, attain a helical configuration in a preferential direction. However, the 
extent of this partial resolution, if it exists, is not known. 


* Contribution No, 1345 from the Sterling Chemistry Laboratory, Yale University. 
+ National Science Foundation predoctoral fellow, 1954-1956. 
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HOMOLOGICAL DIMENSION IN NOETHERIAN RINGS 
By Maurice AUSLANDER AND Davin A. BucHSBAUM* 
DEPARTMENTS OF MATHEMATICS, UNIVERSITY OF MICHIGAN AND PRINCETON UNIVERSITY 


Communicated by A, A, Albert, November 15, 1955 


1. Introduction.—It is our purpose in this note to present some new results in the 
theory of commutative Noetherian rings (with special emphasis on local rings) 
which we have obtained using the concepts and techniques of homological algebra, 
recently introduced by H. Cartan and 8. Eilenberg. Since a more detailed account 
of these and other results will soon appear elsewhere, we take the liberty of omitting 
most proofs and bibliographical material. 

2. Preliminaries.—Throughout this note, R will denote a commutative, No- 
etherian ring. If J is an ideal in R, we denote the rank of J by dim /. We define 
dim R = rad dim 7, where J runs through all the ideals of R. By the integer [/], 


we mean the smallest number of elements which generate 7. We say that R is a 
local ring if R has a unique maximal ideal which we shall always denote by m._ It 
is a classical result that dim R = dim m < [m]. R is called a regular local ring if 
dim R = [m]. 

Let M be an R-module. <A projective resolution of M is a projective, acyclic 
complex X over M, i.e., a sequence of R-modules and R-homomorphisms, 


dn 
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such that each X, is a direct summand of a free R-module, d,.1(X,41) = d,~'(0) for 
alln >0, and Xo/d\(X1) AM. We define the homologieal dimension of WM, h dim 
M, to be the smallest integer n such that there is a projective resolution of M with 
Xn = 0, if such an integer exists. Otherwise we define the homological dimension 
of M tobe ~, 

A basic result concerning h dim M is that if 0 M” is an exact 
sequence of R-modules, with M projective and /” not projective, then h dim M” = 
1+hdim MW’. 

The global dimension of the ring Rk, gl dim RF, is defined to be the sup h dim M, 

M 


where M runs through all R-modules. By a result of Auslander, we know that we 
may confine our attention to finitely generated R-modules. 

3. Codimension.—Let R be a local ring with maximal ideal m. Cartan and 
Filenberg have shown that if M is a finitely generated R-module, then h dim M < 
h dim R/m. Therefore, we have the result that gl dim R = hdim R/m = 1 + 
h dim m. 

3.1. Lemma. Let I be an ideal in the local ring R, and X an element in m such 
that I:(x) = I, where I:(x) = eR/rve Tt. Then h dim UI, x) = 1 + h dim 

Suppose now that F is a regular local ring of dimension n. Then R is an integral 
domain, and if m = (wu, ..., Un), (Wi, ..., isa prime ideal for each i = 1,..., 
n. Therefore, from Lemma 3.1, we obtain 

3.2. THeoreM. J/f R is a regular local ring, then gl dim R = dim R. 

From Lemma 3.1 we can also deduce that m belongs to an ideal / in a regular local 
ring if and only if r dim J = hdim m. 

Let J be a proper ideal in the ring R. We say that 1,...,2, € R isan R/I- 
sequence if [:(x,) = I, UI, = m,...,2,) for all = 1, 

.,p—l,and (U,m,...,2,) # R. If Risa local ring, it readily follows that 
the number of elements in an #R/J-sequence is bounded and that m,..., 2,18 a 
maximal R/J-sequence if and only if m belongs to (7, 7, ...,2,). It is trivial to 
see that every R/J-sequence can be imbedded in a maximal one. 

3.3. THrorem. If R is a local ring and I is an ideal in R, then all maximal R/I- 
sequences have the same number of elements. 

Proofs: Suppose that FR is a regular local ring and dim R = n. Then x,..., 
x, isa maximal ?/J-sequence if and only if m belongs to (J, 71,...,2,), or, equiva- 
lently, h dim (J, 7,...,2,) = hdim m. But, by Lemma 3.1, h dim (U/, m,..., 
ty) = p+hdim/. Thus p = hdim m — h dim J, which proves the invariance 
of p. It is easy to see that the theorem is also true for all factor rings of regular 
local rings. If R is an arbitrary local ring, we know that x, ..., 2, is a maximal 
R/I-sequence if and only if 2, ..., 2, is a maximal R*/J*-sequence, where R* is 
the completion of R, and 7* = R*/. Therefore, it suffices to prove the theorem 
for complete local rings. But every complete local ring is a factor ring of a regular 
localring. Q.E.D. 

An interesting consequence of Theorems 3.2 and 3.3 is that every regular local 
ring of dimension 2 is a unique factorization domain. 

Theorem 3.3 leads us to introduce the following new invariants. We define the 
codimension of R/T (codim R/T), where J is an ideal in the local ring R, to be the 
number of elements in a maximal R//-sequence. We can show that codim R = 
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sup h dim VW, where WV ranges over all finitely generated R-modules of finite homo- 
M 


logical dimension. 

Note: Let R bearing, and WM an R-module. We define an M-sequence to be a 
sequence 2,..., «, of elements in R such that 2x; is not a zero division for M, and 
244, 18S not a zero division for M/(m,...,2,)M, foralli =1,...,p— 1. This 
is a direct generalization of an R/J-sequence. 

4. Regular Local Rings.—Let R be a ring and S a multiplicatively closed subset 
of R not containing 0, and let Rs be the ring of quotients of R with respect to S. 
Since Tor,® (Rs, M) = 0, forall n> 0, we have h dimg M >hdimg, Rs If 
isan ideal in Rs, then AR, @,R/RAT. Thus we have gl dim R > gl dim 
R,. Therefore, if R is a regular local ring and P is a prime ideal in R, then gl dim 
Rp is finite. It is well known that for regular geometric local rings, as well as for 
regular, nonramified complete local rings, Rp is also a regular local ring for every 
prime P. This observation, together with some direct computations, led the 
authors to conjecture 

4.1. THrorem. A local ring R is regular if and only if gl dim R is finite. 

The proof of this theorem is a direct result of the following lemmas, the first of 
which was established by the authors and the second of which is an immediate 
consequence of a result of Serre. 

4.2. Lemma. Jf Risa local ring of finite homological dimensions, then gi dim R = 
dim R. 

4.3. Lemma. If R isa local ring with maximal ideal m, then hh dim R/m > {m]. 

4.4. Corotiary. Jf R is a regular local ring, then so is Rp for every prime ideal 
Pink. 

The last result leads one to say that a ring R is regular if and only if fp is a regular 
local ring for all prime ideals P of R. For a regular ring R, we have gl dim R = 
dim R. If gl dim R is finite, it is clear that R is regular. Also, it is not difficult 
to show that all regular rings are finite direct sums of regular integral domains. 

An ideal J in a regular ring F is said to be perfect if h dim R/J = dim J. Since 
hdim, R/I 2 ed dim P, where P ranges over all the prime ideals belonging to J, a 


perfect ideal is unmixed. We are able to establish the following version of the 
Cohen-MacAulay theorem: 

4.5. Turorem. If Ris aregular ring, and I = (x,...,2,) 1s an ideal of rank 
n, then the x1, ... , £, are an R-sequence regardless of the order in which they are 
written. Consequently, I is a perfect ideal and is therefore unmixed. 

It can also be shown that if'dim R < 2, and R is a regular domain, then every 
minimal prime ideal in # is invertible. 

We conclude this note with 

4.6. THrorem. If R is a Zariski ring and M is an R-module, then h “— M> 


h on M*, where R* is the completion of R and M* the completion of Me. If Ris 
a roel ring, then the above inequality becomes an equality and gl. dim R = gl. dim R*. 


* Work done while one of the authors was a National Science Foundation postdoctoral fellow. 


GENERAL THEORY OF ELLIPTIC EQUATIONS* 
By Leon EnRENPREIS 
INSTITUTE FOR ADVANCED STUDY 


Communicated by M. Morse, November 4, 1955 


Let D be a distribution of compact carrier (support).! We shall study the ques- 
tion of when all distributions S which satisfy the equation 


D,S == (1) 


are real analytic functions (or indefinitely differentiable functions). In case D 
is a partial differential operator with constant coefficients, we shall give necessary 
and sufficient conditions on D in order that all the solutions of equation (1) be 
analytic (or indefinitely differentiable). Partial results in this direction were ob- 
tained by Petrowsky*? and by Hérmander.* In the general case where D is a dis- 
tribution of compact carrier, our results are not quite complete. However, we 
show that if all solutions of equation (1) are analytic, D is the composition of a 
translation with a partial differential operator, and then our conditions on D are 
necessary and sufficient. 

Our method of proof of the sufficiency of the conditions is the construction of an 
elementary solution with the desired properties.! This elementary solution is con- 
structed by means of the Fourier transform. The necessity of our conditions will 
be shown by exhibiting certain exponential sums which define distributions S satis- 
fying equation (1). 

We wish to remark that our method of deriving the sufficiency of our conditions 
by means of the elementary solution has a serious defeet which explains why our 
results are iot complete: there may exist distributions D of compact carrier for 
which all solutions of equation (1) are indefinitely differentiable but such that D 
has no elementary solution.‘ 

Let R denote real Euclidean space of dimension n; C is complex Euclidean n-space. 
By C! we denote the space of all z « C such that 9(z;) 2 Ofor ally. For / any rota- 
tion of C which maps FR into R, we write C' for IC’. Let D be a distribution of 
compact carrier, and denote by P the Fourier transform of D, so that P is an entire 
function of exponential type on C which is slowly increasing on R. Let V denote 
the set of zeros of P in C; assume that V 9 Ris compact. We want to define an 
elementary solution for D in some sense as the inverse Fourier transform of 1/P. 
There are two difficulties in doing this: the compact set of zeros V 9 R and the 
convergence difficulty at infinity. 

1. The sete V © R.—It can be shown that, if A is any open wiaineiliiia of 
V* (a cube in R containing V 9 R) in C, then the map PF +fy« F(z) dz is acon- 
tinuous linear function on PA,, where Ax, is the space of functions analytic on K 
with the topology of uniform convergence on the compact sets of A, and PA,g = 
| PF} max With the topology induced by Ax. Using the Hahn-Banach theorem, 
we can find a measure u on C of compact carrier such that, for F ¢ Ax, {PF du = 
F(z) dz® Suppose that fp—y+{1/P(z)] exp dz converges in the space of 
distributions on R to a distribution 7 (say). Then it is readily verified that the 
sum of 7’ and the inverse Fourier transform of u is an elementary solution for D. 
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Since wu is a measure of compact carrier on C, the inverse Fourier transform of u 
is an entire function of exponential type and thus is as regular as we want. 

2. The set R — V*.—For each x ¢ R, denote by p(x) the distance from x to V. 

Condition a: lim log 

p(x) 

Condition b: There exists ad > 0 such that, for each r > 0 and for each 1, we can 
find a real n-chain -! ¢ iA with the following properties: 

(a) The boundary of - pind te boundary of V*. 

(8) The element of : area on | 4dx,...dx, for | large. 

(y) P-\(z)| = (0( Alar chain in C lying between —! and R. 

(6) For 2 —!, the distance from z to R is 2 r log |z| for | z| large. 

It is easily seen that if P satisfies conditions a and b, then the integral fg—y« 
exp (7z' )[{1/P(z)] dz converges in the space of distributions on FR and is equal to 
S = exp (iz')[1/P(z)] dz for any land Now suppose that x eR, x = (a1,..., 


rn), t) > 0, for allj. Forze —}, if | 2| is large, we must have, for some k, 9(z,) = 
log Thus 
n 


(iz'x) exp (—(z,)a,) S (2) 


From relation (2) we deduce easily that, for any a > 0, we can find an r so large 
that f= exp (ix’){1/P(z)] dz is a times differentiable in a neighborhood of x. 

The case x € I, x ¥ 0, not - x, > 0, may be easily reduced to the above by a 
linear transformation. 

THeoreM |. Let D satisfy conditions a and b. Then all the distribution solutions 
of equation (1) are indefinitely differentiable functions. Conversely, if all the distribu- 
tion solutions of equation (1) are indefinitely differentiable functions, then D satisfies 
condition a. 

Proof: If D satisfies conditions a and b, then the above implies that there exists 
an elementary solution for D which is an indefinitely differentiable function outside 
the origin. The fact that all solutions of equation (1) are indefinitely differentiable 
functions then follows without difficulty.! 

Suppose that D does not satisfy condition a, and let® F ¢ D be so chosen that 
F(0) = 1. We can construct a sequence of points b’ « V such that |b’ + 1) — 
| > 1 for all j, and, for each 1, < 1/2. Morevoer, for some 


jal 
a> Oand any k, < < a log where we have written = (b/,..., b,’). 
We find easily that the series > j exp (ib’") converges in the space 9’ of distribu- 
tions on R, say U = ¥) j exp (ib’*). It is clear that Ds U = 0. On the other 
hand, U cannot belong to the space & because, if f denotes the inverse Fourier 
transform of F, then the set fexp (—ib’-)f} is bounded in &’, but 


U- [exp = — 


Thus U ¢ &, which completes the proof of Theorem 1. 
If D has an elementary solution, we can show by use of the results of a previous 
paper‘ that condition b is a consequence of condition a. We have thus 
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THeoreM 2. Jf D has an elementary solution, then all the solutions of equation (1) 
are indefinitely differentiable functions if and only if D satisfies condition a. 

We shall now study the analytic case. 

Conditionc: For some ¢ > 0, p(x) = ¢| for |x! large. 

Lemma. /f D satisfies condition c, then the carrier of D is one point, i.e., equation 
(1) is the composition of a translation with a partial differential equation. 

Proof: By multiplying P by functions of the form P(2, ... , 2-1, —2), 
z,), We obtain a function P* which is invariant under z, > —z, for any k. By the 
theorem of addition of carriers, it is sufficient to show that P* is a polynomial. 
For fixed 2, , 2j-1, 2j+1, +, Zn, We see that the function of one complex vari- 
able, 2; > P*(a, , 2), Zn), Satisfies condition ¢ and thus must be a 
polynomial.’ The result then follows. 

By reasoning similar to that used in the proof of Theorems | and 2, we can deduce 

THeoreM 3. A necessary and sufficient condition that all distribution solutions of 
equation (1) be analytic functions is that D satisfy condition c. If this is the case, 
then equation (1) is the composition of a translation with a partial differential equation 
with constant coefficients. 

Example: For k = 2, it is easily seen that the polynomial in two letters 2°" + 
2.” satisfies condition a but not condition e. 

The above can be extended without difficulty to systems of equations. More- 
over, the classes ‘analytic functions” and “indefinitely differentiable functions” 
can be replaced by suitable classes of quasi-analytic and non quasi-analytic 
functions. 


Added in proof: 1. In case D does not have an elementary solution, we can still show that, 
if D satisfies condition a, then there exists a linear function LZ on a suitable space of non-quasi- 
analytic functions, such that D * Z = 6, and such that L is an indefinitely differentiable function 
outside the origin. We can then deduce that condition a is necessary and sufficient in order that 
all distribution solutions of (1) be indefinitely differentiable functions. 

2. In case D is a partial differential-difference operator which satisfies condition a, then we 
can show that the support of D is one point, so that equation (1) is equivalent to a partial dif- 
ferential equation with constant coefficients. 

* Work supported by National Science Foundation Grant NSF 5-G1010 at the Institute for 
Advanced Study. The author is on leave from the Johns Hopkins University. The present note 
is asummary of results and methods; the details will appear elsewhere. 

1 Schwartz, Théorie des distributions, Vol. 1 (Paris, 1950). 

2 “Sur l’analyticité des solutions des systémes d’équations différentielles.”’ Mat. Shornick, 5, 
3-68, 1938. 

“On the Theory of General Partial Differential Equations,” Acta Math., 94, 161-248, 1955. 

4See L. Ehrenpreis, ‘Completely Inversible Operators,’’ these PRocEEDINGs, 41, 945-946, 
1955, where a necessary and sufficient condition for the existence of an elementary solution is given. 

5 See L. Ehrenpreis, “Mean Periodic Functions. I,’ Am. J. Math., 77, 293-328, 1955, for the 
properties of the Fourier transform on the dual spaces of spaces of analytic functions. 

® ) is the space of indefinitely differentiable functions of compact support on R; D’ is the space 
of distributions on R; D is the Fourier transform of D; & is the space of indefinitely differentiable 
functions on R; &’ is the dual of &. 

7 See L. Schwartz, “Théorie générale des fonctions moyenne-périodiques,’’ Ann. Math., 48, 880, 
1947. 


SPECTRAL SYNTHESIS FOR THE CANTOR SET* 
By Cari 8. Herz 
DEPARTMENT OF MATHEMATICS, CORNELL UNIVERSITY, ITHACA, NEW YORK 
Communicated by S. Bochner, November 17, 1955 


Let C denote the Banach space of bounded continuous functions on the real line 

with the norm, | g'| = sup | g(x)|. The spectrum! of a function g « C is the 
x 

closed set A(¢) consisting of those points tsuch that feL' and Sele — y) fly) dy = 

0 imply F(t) = 0, where F is the Fourier transform of f, F(s) = (29)~' f exp 

(—isx) f(~) dx. Every closed subset of the real axis is the spectrum of some func- 

tion in 

The spectral synthesis problem is this: When is it true for a closed set 7 that 
geCwith A(g) TandfeL' with F(t) = 0, forte T,imply g(a — y) f(y) dy = 
0? Spectral synthesis on 7 is equivalent to the following statement: If / is a 
closed ideal in the Z' convolution algebra over the reals and 7’ is the set of zeros 
common to the Fourier transforms of functions in /, then every function in L! 
whose Fourier transform vanishes on 7 belongs to /. It is well known that spectral 
synthesis holds if the boundary of 7 contains no perfect set. Our aim is to show 
that the same is true when 7’ is Cantor’s ternary set. The proof is also applicable 
to other Cantor sets, including ones with positive measure. 


Given h > 0, let = sin? he when #0, K,(0) = 1. For 


form 
#,(t) = exp (—itx) Ky(x) g(a) da. 


,(t) = O unless ¢ is at a distance less than h from A(g).? Hence, if A(¢) is compact, 
®, vanishes outside a compact set, and the Poisson summation formula* may be 
applied to yield 


+ + = ®, (nh) exp (inhx). (1) 


The sum on the left converges absolutely and uniformly in every compact set to a 
function g(x); the sum on the right is finite.t In particular, taking ¢ = 1 in 
equation (1), )> Ka(« + 2xh-'m) = 1. From this it is easy to see that, for any 
< and g(x) > g(x) ash—~0. Thus, for any f L', 


Jim S — y) Sy) dy = — y) Sy) dy. (2) 


If ¢ has a compact spectrum and F is the Fourier transform of f, it follows from 
equation (1) that 


en(x — y) fly) dy = h &,(nh) F(nh) exp (inha). (3) 


Now let 7’ be Cantor’s ternary set. We have to show that if A(g) ¢ 7 and 
F(t) = 0 fort « 7, the right-hand side of equation (2) is zero. Let h tend to zero 
through the values 3~? as p— ©. A number of the form n+3~” either belongs to 
T or is at a distance greater than or equal to 3~” from it. Hence, whenever h = 


42 
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3-”, either F(nh) = 0 or ® (nh) = 0. Thus the sum on the right-hand side of 
equation (3) is identically zero. Combining equation (3) with equation (2), we 
have the desired result. 


* The research for this paper was supported by the United States Air Force under Contract 
No. AF18(600-685) monitored by the Office of Scientific Research. 

1 A discussion of the problem along with the necessary references will be found in Harry Pol- 
lard, ““The Harmonic Analysis of Bounded Functions,’’ Duke Math. J., 20, 499-512, 1953. 

2 See thid. 

3 See S. Bochner, Fouriersche Integrale (Leipzig, 1932), p. 33. 

4 The Levitan polynomial. See N.I. Achieser, A pproximationstheorie (Berlin, 1953), p. 146. 
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Introduction.—This paper presents a central limit theorem for a sequence of 
dependent random variables 


The assumptions required are the usual assumptions on second and 2 + 6 order 
moments and a strong mixing condition. The theorem is of interest for two reasons. 
All general central limit theorems for dependent random variables formalize in 
some sense a heuristic notion of A. Markoff to the effect that one expects a central 
limit theorem to hold for X,, Xo, ... , if the random variables behave more like 
independent random variables the farther they are separated (assuming that ap- 
propriate moments exist). An interesting discussion of this intuitive notion is 
given in 8. Bernstein’s paper on the central limit theorem.' The strong mixing 
condition used in this paper seems to be a more intuitively appealing formalization 
of this notion than most others. The condition is also of interest because it is a 
strong version of the mixing condition encountered in ergodic theory (see Hopf,? 
p. 35). 

The Mixing Condition. The mixing condition is assumed to hold for sets of the 
form 


where k; <... < ks. Given any two such sets A and B, introduce the following 
notion of distance d(A, B) between them. Consider the smallest interval contain- 
ing the indices of the random variables in terms of which A is defined. Call this 
interval /(A). In the case of set (1), this interval is [h;, k,]. The distance d(A, B) 
is set equal to the distance between the intervals /(A) and/(B). Note that d(A, B) 
does not satisfy the usual properties of a distance; for example, one may have d(A, 
B) > d(A, C) + d(B, C) for some set C. The strong mixing condition is satisfied 
by the sequence of random variables X,, Xo, ... , if 
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|P(A B) — P(A) P(B)| < B)) (2) 


for all sets A, B of form (1), where f(n) is a function defined for n = 0, 1, 2,.... 


and decreasing to zero asn—> ©. 
The Central Limit Theorem..Assume that the mean values LX, = 0, as one 
could in any case subtract the mean values HX, from the random variables X,,, 


respectively. The assumptions are 


— 0) 
=a 


as b — a— ~, where h(m) 7 © asm—> 


= o(h(b — a)! +6? 


as b — a— o for some 6 > 0. 
2. The strong mixing condition (2) is satisfied by the sequence of random 


variables X,, Xo, ... 
Introduce the following notation: 


= Xj, 


j=1 
rpn + (r — 


X;, 
=(r — 1)(tn + aq) +1 


r(pn_+ Qn) 
+ 
where k(p, + dn) = n. The numbers k, pp, qn Will be chosen in such a way that 
Pns Qn; k—> and O0asn—> 


We expect to show that 


is asymptotically normally distributed. Now 
| > V 
| 
h(pn))” 


Choose the sequences gn, Pn», k such that 


0. 
h(pn) 


k V, 


V kh(pn) 


4 
A 
Ei X; 
Ss 
= 1 
Then 
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in probability as n + ©. It can therefore be neglected in further computations. 
We shall later see that this can be done. 
Set 


and let the set 


U, 


where m, is an integer. Note that 


(m +...+ mk +k) r= 


Vich(p,) 
A(r, n, m,, 


Lemma 


if te = 
Now 


The result is obtained by an application of the Tchbycheff inequality. 
LEMMA 2 


| k 
| A At, n, M,, Il P(A(r,n, m,, 


+...4+ < x (mi +...+ m)icsxr=1 


+ 


k 
The probability contributed by all the sets f¥ A(r, n, m,, 6) for which max 
rail 
k 
| U,/ V ich(p,)| > 4 is at most «. Consider the sets f} A(r, n, m,, 6) for which 


max | U,/+/ kh(pr)| < te. By repeated application of condition (2), one can see 
that 


k k | 
P ( NM Alr, i) — II P(A(r, n, m,, < kf (qn): 
| r=1 
Since there are (2t,/5)* sets of this form, we obtain the desired inequality. 
Let 


* Ge, (2) 


G,, n(x) < 
r, nt) = 
Vih(p,) 
i 
(my +... 4+ 
U, 
max >t) < 
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be the convolution of G,, ..., Note that 


k 
* Gy, n(X) < II P(A(r, n, M,, 5)) 
(m +...+ meisx (4) 


+8) 
and 


k 
Gy,» Gy, n — bb) < Il P(A(r, n, m,, 5)) 
(m+...+ m+k)6 5x r=1 (5) 
Now , * .. . * G, ,(x) tends to the normal distribution as n © if k—~ o, 
We want to let hk, pn, Gn > ©, kpn ~ n, and ké > 0 in such a way that 


‘ k 
k > 0 (6) 


KR(Qn) 
h(pn) 
Set 6 = 1/k2. Thenké—~0. Since 


0. (7) 


it follows that 


‘ k k 
€ 


where C = 2/e*. The function f can always be chosen in such a way that f(n) > 
1/nforalln. Ifkis chosen so that 


k < [—log f(qn)]", (8) 


clearly condition (6) will be satisfied. 
Now 


n 
h(n)~E < (‘) 
by the Minkowski inequality. But then 

~2 
Mg 
h(n) ~ 


h(n) ~ 
Condition (4) will then be satisfied if 


h 
h(qn) = 0 (9) 
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Sequences k, p,, Gn, > ©, kp, ~ n can be chosen so that conditions (8) and (9) are 
satisfied and hence conditions (6) and (7) are satisfied. On using inequalities (3). 
(4), and (5) and Lemma 2, we see that 
k 
U, 


is asymptotically normally distributed. Since 
k 
V, 


— 0 


in probability as n > ©, it follows that 
S, 


is asymptotically normally distributed. 

Remarks.—The result obtained includes the result of Héffding and Robbins.* 
It would be of very great interest to see how much stronger the notion of a strong 
mixing condition is than that of an ordinary mixing condition in the case of a 
strictly stationary process. 


1S. Bernstein, “Sur l’extension du théoreme limite du calcul des probabilités aux sommes de 
quantités dépendentes,” Math. Ann., 97, 1-59, 1927. 

2 I. Hopf, Ergodentheorie (““Ergebnisse”’ series). J. Springer, Berlin, 1937. 

3 W. Hoffding and H. Robbins, “The Central Limit Theorem for Dependent Random Vari- 
ables,’’ Duke Math. J., 15, 773-780, 1948. 
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